Yang Luo,??** Lu Lu,?** Alexandra H. Wortley,* 
De-Zhu Li,?** Hong Wang,?** and Stephen 
Blackmore* 


EVOLUTION OF ANGIOSPERM 
POLLEN. 3. MONOCOTS' 


ABSTRACT 


Monocots, with ca. 65,000 species in 78 families and 12 orders as classified by the Angiosperm Phylogeny Group 
(Angiosperm Phylogeny Group IL, 2009), have traditionally been known as a distinct assemblage from dicots within 
angiosperms and are now seen as a lineage that diverged from within the more basal of the angiosperm groupings. In this study, 
the third in a series analyzing pollen characters across the angiosperms, we illustrate the pollen morphological diversity of 
monocots and analyze 19 palynological and two ecological characters for 120 taxa in 71 families covering the 12 monocot orders 
and 16 taxa of nine orders in basal angiosperms. Pollen morphological data from previous works and our investigations were 
optimized onto a new maximum likelihood tree reconstructed from an existing DNA matrix of Chase et al. (2006) using Fitch 
parsimony, maximum likelihood, and hierarchical Bayesian analysis. From these analyses we infer evolutionary patterns in 
palynological characters, assess their systematic value, and investigate two aspects (pollination type and habitat moisture) of 
their ecological adaptation. The highest levels of pollen variation were shown to exist in the Alismatales and Commelinales, with 
lower levels seen in the Asparagales, Dasypogonales, and Zingiberales; the most variable characters across the phylogeny were 
found to be pollen outline in polar view, size, and tectum extent. We infer unambiguous plesiomorphic states for monocots and 
report significant transitions in character states at various levels within the monocot assemblage. Analyses of correlated 
evolution reconfirmed the hypothesis of an association between exine reduction (or complete loss) and habitat moisture and 
found significant correlations between various states of exine loss and habitat for hydrophytic-helophytic plants. The presence or 
absence of the exine itself was found to be more significantly correlated than individual structures, in terms of association with 
hydrophytic—helophytic plants. The most rapid rate of state changes in pollen characters, in the evolutionary history of monocots, 
is estimated to have occurred during the Albian-Turonian stage; our work may provide insights into the identification of 
enigmatic fossil pollen grains from this geological time. 


Key words: Character evolution, correlated evolution, exine stratification, habitat moisture, monocots, pollen morphology, 
pollination type, systematic significance. 


Monocotyledonous plants have been recognized as 
a distinct group, based principally on the single 
cotyledon of the seed, since the 18th century (Ray, 
1703). Monocots are one of the most distinctive major 
lineages of angiosperms, with ca. 65,000 known 
species in 78 families and 12 orders, representing ca. 
25% of the total diversity of angiosperms (Angio- 
sperm Phylogeny Group III, 2009; Reveal & Chase, 
2011). Cladistic studies based on morphological data 
have indicated putative synapomorphies for monocots 


that include an adventitious root system, parallel- 


veined leaves, scattered vascular bundles in the 
shoots, cambium presence, and successive micro- 
sporogensis (Donoghue & Doyle, 1989a, 1989b; 
Loconte & Stevenson, 1991; Doyle & Donoghue, 
1992; Doyle & Endress, 2000). Tomlinson (1995) 
stated that the uniqueness of the monocot vascular 
system among seed plants and the lack of homology 
in organization between the primary vascular systems 
of monocots and other angiosperms might suggest an 
independent origin of monocots from other angio- 


sperms. However, Tomlinson's observations are 
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inconsistent with other morphological characters and 
molecular data. Scattered vascular bundles and 
monocotyledonoid floral features are found elsewhere 
in angiosperms, among basal lineages (Soltis et al., 
2005), and the monocots are now known to nest 
within angiosperms, closely related to the magnoliids 
(a clade comprising Canellales, Laurales, Magno- 
liales, and Piperales), as supported by analyses of 
several nucleotide sequences (Chase et al., 1993; 
Soltis et al, 1997; Qiu et al., 1999, 2000, 2005, 2006, 
2010). Thus, monocots currently form one of the best 
supported major clades within angiosperms (Chase, 
2004; Chase et al., 2006). 

Phylogenetic studies of the monocots have made 
immense strides over the past 20 years, through 
increasingly extensive taxon sampling and the 
application of molecular data from plastid, nuclear, 
and mitochondrial genomic regions (e.g., Chase et al., 
1993, 2006; Duvall et al 1993a, 1993b, 2006; 
Bremer, 2000; Davis et al., 2004; Tamura et al., 
2004; Givnish et al., 2006, 2010; Graham et al., 
2006; Petersen et al., 2006; Li & Zhou, 2007). The 
most basally branching monocots have been consid- 
ered, from morphological inference, to be the aquatic 
alismatids (a clade including all families of Alisma- 
tales except Araceae and Tofieldaceae; Hutchinson, 
1934; Cronquist, 1968, 1981; Takhtajan, 1969, 
1991; Stebbins, 1974; Thorne, 1976), Dioscoreales 
(Dahlgren et al., 1985) or Melanthiales (Thorne, 
1992a, 1992b). Since some of the first molecular 
studies of monocots (Chase et al., 1993; Duvall et al., 
1993a, 1993b), Acorus L. has been inferred as the 
basalmost branching taxon and sister to all other 
monocots, a finding that has received support from 
almost all subsequent DNA-based analyses. 

The most recent angiosperm classification provid- 
ed by the Angiosperm Phylogeny Group III (2009) 
recognized 12 monocot orders: Acorales, Alismatales, 
Asparagales, Dioscoreales, Liliales, Pandanales, 
Petrosaviales, and the commelinid orders Arecales, 
Commelinales, Dasypogonales (Dasypogonaceae, An- 
giosperm Phylogeny Group III, 2009, here recognized 
as an order following Reveal & Chase, 2011), Poales, 
and Zingiberales. In this scheme, Acorales are the 
first branch to diverge from the rest of the group, with 
the Alismatales second, followed by Petrosaviales 
and Dioscoreales-Pandanales. Liliales are sister to a 
large clade of all other core monocots (i.e., 
Asparagales and the commelinid orders). Although 
this phylogenetic framework for monocots is mostly 
well supported at the ordinal level by nuclear, 
plastid, and mitochondrial data (Chase et al., 2006), 
the relationships among some lineages have remained 
controversial. For instance, the positions of the 
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Dasypogonales and Arecales, relative to other 
commelinids, are ambiguous, and relationships 
remain unresolved in Zingiberales, in particular the 
placements of the Heliconiaceae and Musaceae, as 
well as in the most basal branches of the Poales (e.g., 
Kress et al., 2001; Davis et al., 2004; Tamura et al., 
2004; Givnish et al., 2006, 2010; Graham et al., 
2006; Petersen et al., 2006; Barrett et al., 2014). 
Pollen can play a key role in studies of angiosperm 
evolution due to its morphological variability and its 
potential to fossilize (Jones & Rowe, 1999). Multiple 
studies have indicated the taxonomic significance of 
pollen morphology in monocots, including investiga- 
tions in Alismataceae (Argue, 1974, 1976), Araceae 
(Grayum, 1992; Weber et al., 1990; Hesse, 2006), 
Arecaceae (Thanikaimoni, 1970; Sowunmi, 1972; 
Ferguson, 1986; Ferguson & Harley, 1993; Harley & 
Baker, 2001), Cyperaceae (Van Wichelen et al., 
1999; Simpson et al., 2003), Dioscoreales (Caddick 
et al., 1998), Eriocaulaceae (Thanikaimoni, 1965; de 
Borges et al., 2009), Liliaceae (Nair & Sharma, 1965; 
Kosenko, 1999), Smilacaceae (Chen et al., 2006), 
and Tecophilaeaceae (Simpson, 1985b). Similarities 
in pollen wall structure provided evidence for a close 
relationship between Lemnaceae and Araceae (the 
former treated as a subfamily within the latter in the 
Angiosperm Phylogeny Group lll; Hesse, 2006), 
supporting the recent molecular results of Barabé et 
al. (2002) and Cabrera et al. (2003). Distinctive 
aperture patterns and tectum sculptures have sup- 
ported the combination of Heterosmilax Kunth and 
Smilax L. 
Rhipogonum J. R. Forst. & G. Forst. from Smilaceae 
(Chen et al., 2006). Pollen characters, such as 


aperture pattern and spine type, are also consistent 


in the Smilaceae and the removal of 


with treatment of the Eriocaulaceae at various 
taxonomic levels (de Borges et al., 2009). Pollen 
evolution across the monocols was previously sum- 
marized in detail by Zavada (1983), who proposed 
evolutionary trends in pollen aperture and wall 
structure parallel to those proposed for other 
angiosperms by Walker (1974, 1976): this suggested 
similar selective pressures on pollen characters in 
monocots and other angiosperms. Increases in the 
number and number of types of apertures in the 
pollen grains of the Alismatales, Arecidae, and 
Liliales, and a reduction of the suleus among many 
monocot groups, are obvious trends that, however, 
have not been investigated in a modern phylogenetic 
context but, so far, only under traditional classifica- 
lion systems such as that of Cronquist (1981). 

The monocots are commonly considered to have 
originated by the Early Cretaceous (146—100 million 
years ago [Ma]) and diversified into the major extant 
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lineages during the Late Cretaceous (100—66 Ma), 
contemporaneously with other basal angiosperm 
groups (Herendeen & Crane, 1995; Bremer, 2000; 
Gandolfo et al., 2000; Chaw et al., 2004; Janssen & 
Bremer, 2004; Bell et al., 2010; Smith et al., 2010). 
Pollen grains of monocots from the Early Cretaceous, 
such as Liliacidites Couper and other types with a 
distinct graded reticulum, have long been regarded as 
evidence of early monocots (Doyle, 1973; Walker & 
Walker, 1984, 1986). On the other hand, Gandolfo et 
al. (2000) questioned all Early Cretaceous records of 
monocots and stated that the earliest unequivocal 
monocot fossils are from the Turonian. The appar- 
ently sparse fossil record of monocots is probably due 
lo this group's fundamentally herbaceous habit and 
the fact that some monocot pollen is hard to 
distinguish from the monosulcate pollen of some 
basal angiosperms (e.g., Magnoliales; Friis et al., 
2006). Friis et al. (2004) argued that the earliest 
proposed unequivocal monocot fossil was Mayoa 
portugallica Friis, Pedersen & Crane, with poly- 
plicate and inaperturate pollen grains, which was 
placed in the Araceae (Pothoideae). However, 
Hofmann and Zetter (2010) suggested that this 
placement may be incorrect, based on similarities 
to the far older Triassic spore taxon Lagenella 
martinit (Leschik) Klaus. The incompleteness and 
insufficiency of the fossil record have limited our 
understanding of the evolutionary links between 
extinct and extant taxa and present obstacles to 
identifying early fossil monocot pollen grains. 
Therefore, determining plesiomorphic pollen states 
and pollen evolutionary patterns from extant taxa in a 
phylogenetic context may provide new insights into 
the palynological characters of the earliest monocots. 

Evidence from phylogenetics has indicated that 
pollinators played an important role in angiosperm 
diversification (Johnson, 2006; Armbruster & Much- 
hala, 2009; Kay & Sargent, 2009). Pollen grains of 
angiosperms may be transported to the stigma by a 
variety of biotic and abiotic mechanisms for which 
pollen structures are adapted (Ackerman, 2000; 
Dobson & Bergstrom, 2000; Lunau, 2000; cf. Edlund 
et al. 2004). Relationships have been observed 
between pollination strategies and pollen features 
(e.g., Skvarla et al., 1978; Grayum, 1986; Linder, 
1998, 2000; Hesse, 2000; Tanaka et al., 2004; 
Sannier et al., 2009; Wang et al., 2014); pollen wall 
stratification may be affected by pollination syndrome 
(reviewed in Hesse, 2000), and reduced exine 
structures may be associated with hypohydrophily 
De, underwater cross pollination, Tanaka et al., 
2004). Similarly, characteristics such as small grain 
size appear to be linked to pollen delivery by abiotic 
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vectors (wind or water; Ackerman, 2000). Aspects of 
the habitat features (e.g., moisture) have also likely 
played a substantial role in the evolution of pollen 
morphology (Ackerman, 2000; Franchi et al., 2002). 
For instance, “harmomegathy,” a term introduced by 
Wodehouse (1935) to describe alterations in form that 
accompany changes in the degree of hydration of 
pollen grains, allows grains to respond and adapt to 
their environment. Both aperture closure and pollen- 
wall folding act to accommodate a decrease in cellular 
volume due to water loss (Katifori et al., 2010). 

The loss or reduction of the exine has been a 
particularly prevalent and interesting feature of 
monocols. Pollen grains with a partially or fully 
reduced exine are typically found in moist environ- 
ments (such as marshes, swamps, and bogs), where 
the pollen can remain viable without investing in a 
protective exine, and especially in the evolution of 
aquatic groups (Wodehouse, 1935). Exineless or 
exine-reduced pollen has been found in a large 
number of monocot families including Araceae 
(Weber et al., 1998; Hesse et al., 1999; Weber et 
al., 1999), Asparagaceae and Cannaceae (Erdtman, 
1952, 1963; Skvarla & Rowley, 1970; Rowley & 
Skvarla, 1986), Costaceae (Stone et al., 1981), 
Cymodoceaceae (Pettitt, 1976, 1981, 1984; Ducker 
et al., 1978; McConchie et al., 1982; Harris et al., 
1994), Haemodoraceae (Simpson, 1989), Heliconia- 
ceae (Kress et al., 1978; Stone et al., 1979; Kress & 
Stone, 1983), Hydrocharitaceae (Pettitt, 1981; Tana- 
ka et al., 2004), Lowiaceae (Kamelina & Raven, 
1997), Marantaceae (Erdtman, 1952; Sharma, 1968; 
Furness & Rudall, 1999), Melanthiaceae (Takahashi, 
1987), Musaceae (Erdtman, 1952; Sharma, 1968), 
Strelitziaceae (Hesse & Waha, 1983; Kronestedt- 
Robards & Rowley, 1989; Rowley et al., 1997), 
Zingiberaceae (Theilade & Theilade, 1996), and 
Zosteraceae (Pettitt & Jermy, 1975; Pettitt, 1976). 
Elsewhere in angiosperms, exineless or exine- 
reduced pollen is almost entirely restricted to basal- 
branching groups such as Annonaceae (Waha, 1987), 
Ceratophyllaceae (Takahashi, 1995), Degeneriaceae 
and Gomortegaceae (Rowley & Vasanthy, 1980; 
Kubitzki, 1981; Hesse & Kubitzki, 1983; Rowley & 
Vasanthy, 1993), Hernandiaceae (Rowley & Vasan- 
thy, 1980; Kubitzki, 1981; Hesse & Kubitzki, 1983), 
and Lauraceae (Rowley & Vasanthy, 1980; Kubitzki, 
1981; Hesse & Kubitzki, 1983; Gabarayeva et al., 
2010), with an exception being the obligate aquatic 
species of Callitriche L. in Plantaginaceae (Osborn et 
al., 2001). For a fuller listing of exineless taxa, see 
Furness and Rudall (1999), and for reviews of aquatic 
laxa in particular, see Cox (1988), Osborn and 
Philbrick (1994), Philbrick and Les (1996), and 
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Cooper et al. (2000). Hesse and Kubitzki (1983) 
claimed that truly exineless pollen is found only in 
the seagrass taxa of Alismatales, Amphibolis C. 
Agardh (Cymodoceaceae), Halophila Thouars (Hy- 
drocharitaceae), and Thalassodendron Hartog (Cym- 
odoceaceae), with other taxa (e.g., Canna L., 
Heliconia L., and Strelitzia Aiton in Zingiberales) 
retaining at least vestiges of exine. However, Cooper 
et al. (2000) stated that fully exineless pollen can also 
be found in some aquatic species of Callitriche. 
Exine reduction has been described (Furness & 
Rudall, 1999) as analogous to aperture formation. Since 
apertures are defined as an opening or thinning within 
the exine, exineless pollen grains may be functionally 
either monoaperturate or omniaperturate (in this case, 
the single aperture covering the entire surface of the 
grain). However, exineless taxa exhibit a wide range of 
developmental pathways, in contrast to the single 
pathway commonly described for aperture formation, by 
the deposition of apertural shields of endoplasmic 
reticulum that block the deposition of sporopollenin. 
Since exine may be laid down on a fibrillar layer 
(primexine) within the callose wall during the tetrad 
period of microsporogenesis, or mediated earlier by the 
plasma membrane and glycocalyx (Takahashi, 1987), 
exineless pollen grains may be generated either by a 
lack of primexine (perhaps as in Ceratophyllum L., 
Ceratophyllaceae, Ceratophyllales [Takahashi, 1995] 
and partially in Canna, Cannaceae, Zingiberales 
[Skvarla & Rowley, 1970; Rowley & Skvarla, 1986)), 
callose (perhaps the case in Callitriche; Osborn et al., 
2001), or exine, or by subsequent loss of exine (e.g., 
Zingiber Mill., Zingiberaceae, Zingiberales; Theilade & 
Theilade, 1996). However, Canna does have primexine 
in the shape of the spines, which do not become exine. 
Also, within the primexine of typical pollen grains, 
sporopollenin only polymerizes at specific sites—the 
parts of the primexine that do not accumulate 
sporopollenin are later dispersed—appearing as empty 
voids in the mature pollen wall, such as the spaces 
between columellae. Frequently, exine loss is accom- 
panied by some degree of compensation in other layers, 
such as the retention of the primexine or the thickening 
and elaboration of the intine (Kress & Stone, 1982; 
Theilade & Theilade, 1996; Franchi et al., 2002). 
The extensive accumulation of pollen morpholog- 
ical data and an increasingly robust phylogenetic 
framework for angiosperms based on advanced 
molecular analyses together allow for more reliable 
reconstruction of morphological character states for 
extant flowering plants and their evolution in a 
phylogenetic context (Endress & Doyle, 2009; 
Wortley et al., 2015). The pollen morphology and 
evolution of monocots are integral to our understand- 
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ing of morphological evolution in the angiosperms as 
a whole. In this paper, we survey variation in pollen 
characters and infer key character state transitions in 
various lineages based on the most up-to-date 
phylogenetic framework, using Fitch parsimony 
(FP), maximum likelihood (ML), and hierarchical 
Bayesian inference (HB), in order to reveal pollen 
evolutionary patterns and trends throughout the 
monocots and to assess the systematic value of pollen 
characters as evidence for hitherto uncertain rela- 
tionships within this group. We also perform 
correlated evolution analyses between pollen traits 
and ecological characters, to provide further insights 
into adaptive evolutionary mechanisms and driving 
forces for pollen character state changes in monocots, 
as corroborated with data on divergence times for the 
lineage from previous works. The integration of pollen 
data and analysis in this study may also facilitate the 
interpretation of early fossil monocot pollen grains. 


MATERIALS AND METHODS 
POLLEN DATA COLLECTION 


Pollen grains from 20 species in 20 families 
covering all 12 orders within monocots were 
examined by LM and SEM to demonstrate the pollen 
variability of monocots. Pollen samples were collect- 
ed from the living collections of the Kunming Botanic 
Garden (KBG) and Xishuangbanna Tropical Botan- 
ical Garden (XTBG), China, and from herbarium 
specimens of the Kunming Institute of Botany (KUN), 
China, Royal Botanic Garden Edinburgh (E), U.K., 
and Australian. National. Herbarium. (CANB; for 
voucher details see Appendix 1). Anthers from 
herbarium material for LM and SEM observation 
were softened in warm water for 5 to 10 minutes and 
then dehydrated through an ethanol series to 10096 
ethanol. The acetolysis method was inapplicable 
because most pollen grains of monocots are fragile 
and liable to collapse. For SEM, anthers were 
mounted on specimen stubs with double-sided sticky 
tape and then carefully dissected to expose the pollen 
grains. The stubs were sputter-coated with gold. 
Pollen morphology was examined using a Hitachi 
(Tokyo, Japan) S-4800 SEM at 10.0 kV (KUN). To 
dry extremely delicate pollen grains (such as those 
with a very thin pollen wall) while preserving their 
surface structure for SEM observation, treatment with 
2,2-dimethoxypropane (DMP) and critical-point dry- 
ing (SPI-DRY Critical Point Dryer, No. 13200-AB 
Manual CPD; SPI Supplies, West Chester, Pennsyl- 
vania, U.S.A.) was performed according to Halbritter 
(1998). Pollen size (the maximum diameter of a single 
grain) was measured under LM. Note that, due to our 
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preparation methods, sizes may not be directly 
comparable with those in the literature derived from 
acetolyzed grains. 


POLLEN CHARACTERS AND CODING STRATEGY 


Pollen data were extracted from the literature 
(primarily SEM and transmission electron microscopy 
[TEM] observations), online databases, and new 
observations presented in this paper (see Appendix 
1). Palynological terminology follows that of Punt et 
al. (2007) and Wortley et al. (2015). 

The character selection and coding strategy 
followed those of Wortley et al. (2015) and Harley 
and Zavada (2000). For reconstruction of character 
evolution, 19 pollen characters were selected for 
coding: dispersal unit, polarity, symmetry (in polar 
view), basic shape, shape class, outline in polar view 
(amb), size (diameter of largest axis), aperture 
number, aperture position, aperture membrane, 
aperture shape, annulus, exine extent, tectum extent, 
supratectal element presence/absence, supratectal 
element shape, tectum sculpture, infratectum struc- 
ture, and foot layer presence/absence (see Table 1 for 
coded character states). Two ecological characters, 
pollination type and habitat moisture, were also 
coded (see Table 1 for coded character states). 
Characters were treated as unordered binary or 
multistate characters, with inapplicable states coded 
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as and unknown data as “?.” We adopted two 
coding methods, the comprehensive method sensu 
Wortley et al. (2015), which includes polymorphic 
states, and the democratic method (coding for the 
most common state; Bininda-Emonds et al., 1998). 
We adopted the composite coding strategy (incorpo- 
rating inapplicable states into a more inclusive 
multistate character) advocated by Maddison (1993) 
for dealing with the inapplicable states. Data matrices 
are presented in Appendices 2 and 3. 

For studies of correlated evolution, seven pairs of 
palynological character states and one pair of 
ecological character states, selected on the basis of 
previous hypotheses of correlation (Kress & Stone, 
1982; Theilade & Theilade, 1996; Ackerman, 2000; 
Franchi et al., 2002; Katifori et al., 2010), were 
coded as binary presence/absence characters (see 
Part II of Table 1 for coded character states). 
Correlations were tested between six exine stratifica- 
tion characters: exine, tectum, eutectum (whether the 
tectum provides an effective seal against water loss), 
infratectum (presence of alveolar, granular, columel- 
late, or structureless material beneath the tectum and 
above the foot layer or endexine), columellate 
infratectum, and granular infratectum, and two 
ecological characteristics: helo-hydrophytes or hy- 
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drophytes (hydrophytes are defined as plants adapted 
to aquatic environments, and helophytes to wet but 


marshy habitats; McDonald et al., 2002). This data 


matrix is presented in Appendix 4. 


ANALYSES OF CHARACTER EVOLUTION 


Of recent phylogenetic investigations focused on 
monocots (Davis et al., 2004; Tamura et al., 2004; 
Chase et al., 2006; Givnish et al., 2006, 2010; 
Graham et al., 2006; Petersen et al., 2006), the 
seven-gene phylogenetic tree of Chase et al. (2006) 
provides the most powerful study of relationships 
across the group to date (Givnish et al., 2010), with 
the broadest taxonomic coverage across the widest 
range of genetic loci (nuclear, mitochondrial, and 
chloroplast regions). The results were highly congru- 
ent with those of other recent phylogenetic studies. 

To enable comparative studies based on a 
phylogenetic tree with branch lengths, a tree was 
generated based on the original DNA matrix of Chase 
et al. (2006), obtained from Mark Chase, containing 
120 taxa in 71 families covering all 12 orders of 
monocots (following the Angiosperm Phylogeny 
Group III system) plus 16 outgroups of nine orders 
of the basal angiosperms, using maximum likelihood 
(ML) inference in RAxML (Stamatakis, 2006), 
including a simultaneous bootstrap analysis of 2000 
replicates, under the GTRGAMMA model. 

Two strategies for optimization of pollen characters 
were used: (1) analysis of the comprehensive dataset 
with Fitch parsimony and hierarchical Bayesian 
inference (abbreviated as CFP and CHB, respective- 
ly) and (2) analysis of the democratic dataset with 
Fitch parsimony, maximum likelihood, and hierar- 
chical Bayesian inference (hereafter abbreviated as 
DFP, DML, and DHB, respectively). The ML method 
as implemented in Mesquite does not accept 
polymorphisms, which were numerous in our com- 
prehensive dataset; therefore, only the democratic 
matrix was analyzed under ML. Ecological characters 
were analyzed only under CFP and CHB. A series of 
CFP analyses were conducted to determine differ- 
ences in pollen evolutionary pattern on the much- 
debated alternative topologies within commelinids 
and Poales (Janssen & Bremer, 2004; Chase et al., 
2006; Graham et al., 2006; Petersen et al., 2006; 
Givnish et al., 2006, 2010; Soltis et al., 2011; Barrett 
et al., 2012; Davis et al., 2012; see Fig. 1), with the 
same number of terminals as that of Chase et al. 
(2006) but in different phylogenetic skeletons as 
above. 

FP and ML (using the Mk-1 model) optimizations 
were performed in Mesquite 2.75 (Maddison & 
Maddison, 2011), and HB optimizations using 
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Table 1. 
study. Part I indicates character states for 19 pollen characters (1-19) and two ecological characters (20 and 21), which are mapped 


Pollen and ecological characters and their states used in analyses of character optimization and as defined in this 


in Figure 7A, B. Part II indicates binary characters of pollen morphology (al to a6) and ecology (bl and b2) that were used in 
analyses of correlated evolution. Part II is not mapped on Figure 7A, B. 


Part I 
1. Dispersal unit: 0, monad; 1, permanent tetrad (includes pseudomonads) 

. Polarity: O, apolar; 1, heteropolar; 2, isopolar (or subisopolar) 

. Symmetry (in polar view): O, bilateral; 1, radial; 2, asymmetrical 

. Basic shape: 0, boat-shaped; 1, globose; 2, filiform 

. Shape class: 0, oblate; 1, spheroidal; 2, prolate 

. Outline in polar view (amb), states as defined by Walker & Doyle (1975): 0, circular (longest equatorial axis [LA |/shortest 
equatorial axis [SEA] = 1); 1, elliptic (1 < LA/SEA < 1.5); 2, oblong (1.5 € LA/SEA < 2); 3, elongate (LA/SEA > 2) 

7. Size (diameter of longest axis), states as defined by Walker & Doyle (1975): 0, small (10-24 pm); 1, medium (25-49 um); 

2, large (50-99 um); 3, very large (100-199 um); 4, gigantic (> 200 um) 
8. Aperture number: 0, zero; 1, one; 2, two; 3, three; 4, more than three 
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9. Aperture position: 0, polar: distal; 1, polar: proximal; 2, equatorial; 3, global 
10. Aperture membrane: 0, smooth (without ectexine); 1, sculptured (with numerous ectexine elements but not operculate); 2, 
with operculum 
11. Aperture shape: 0, colpate (when at the distal or proximal pole of a pollen grain, we usually term it sulculate); 1, porate; 2, 
zonate; 3, spiral; 4, syncolpate 
12. Annulus: 0, absent; 1, present 
13. Exine extent: 0, none or merely small patches appearing as sculpture (exineless); 1, homogeneous exine (atectate); 2, fully 
present 
14. Tectum extent: 0, intectate (without a tectum but with sculpturing); 1, semitectate; 2, eutectate 
15. Supratectal elements: 0, absent; 1, present 
16. Supratectal element shape: 0, gemmate; 1, echinate; 2, verrucate 
17. Tectum sculpture: 0, imperforate; 1, perforate (micro- and macro-); 2, reticulate (foveolate); 3, rugulate; 4, striate; 5, 
fossulate; 6, areolate; 7, cupulate 
18. Infratectum structure: 0, columellate; 1, granulate; 2, intermediate between granulate and columellate. Terminology is taken 
from Doyle (2005). 
19. Foot layer: 0, absent; 1, present 
20. Pollination type: 0, self-pollination; 1, wind pollination; 2, insect pollination; 3, water pollination: 4, bat pollination; 5, bird 
pollination; 6, non-flying mammal pollination 
21. Habitat moisture; states as defined by McDonald et al. (2002): 0, xerophytic (adapted to arid and low water availability 
soils); 1, mesophytic (adapted to soils of moderate moisture content); 2, helophytic (adapted to wet but marshy habitats); 3, 
hydrophytic (adapted to aquatic environments). 
Part II 
al. Exine: 0, absent; 1, present 
a2. Tectum: 0, absent; 1, present 
a3. Eutectum: 0, absent; 1, present 
a4. Infratectum: 0, absent; 1, present 
a5. Columellae: 0, absent (all other infratectum types); 1, present 
a6. Granules: 0, absent (all other infratectum types); 1, present 
bl. Helo-hydrophytes: 0, absent (all other habitat types); 1, present (plant adapted to wet and marsh or aquatic habitat) 
b2. Hydrophytes: 0, absent (all other plant types); 1, present (plant adapted to aquatic habitat) 


BayesTraits 1.0 (available from ` «http://www. In the discussion that follows, the state of the most 


evolution.rdg.ac.uk/BayesTraits.html>). Both CHB 
and DHB analyses were run in BayesTraits using 
BayesMultistate with empirically selected, most 
appropriate settings: a reversible-jump hyperprior 
model, prior distribution range for uniform seeding 
distribution (0—10 distribution, 0-20 distribution, or 
0-30 distribution depending on the character), rate 
variations ranging from 5 to 35, a sampling frequency 
of 300 to 800 generations, a burn-in of 10,000 
generations, and a total of 5,000,000 generations. 


recent common ancestor (MRCA) was presented only 
for derived states, in contrast to the plesiomorphic 
states of the monocots. For ML and HB analyses, we 
took the state with the highest likelihood or 
probability at each node to represent the state of this 
node on the phylogeny. 


TESTS OF CORRELATED EVOLUTION 


BayesDiscrete tests for correlated evolution on 12 


pairs of character states from the palynological and 
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branching group of Poales). A indicates the Arecales; D, Dasypogonales; P, Poales; CZ, Commelinales—Zingiberales; Bro, 
Bromeliaceae; Typ, Typhaceae; and RCXRG indicates Rapateaceae and the clade comprising the cyperids, xyrids, restiids, and 
graminids within Poales. The nine topologies correspond to the following recent phylogenetic studies: —A. Janssen & Bremer 
(2004, p. 390), maximum likelihood (ML) tree in this study. —B. Chase et al. (2006, seven-gene maximum parsimony [MP] tree, 
pp. 67-68), Peterson et al. (2006, four-gene summarized tree, p. 56, fig. 2A). —C. Chase et al. (2006, four-gene MP tree, pp. 
69—70). —D. Givnish et al. (2010, ML tree, p. 592), Barrett et al. (2012, pp. 9, 12), Davis et al. (2012, p. 329). —E. Graham et 
al. (2006, p. 6), Givnish et al. (2010, MP tree, p. 590). —F. Givnish et al. (2006, p. 39). —G. Chase et al. (2006, seven-gene MP 
tree, pp. 67—69), Givnish et al. (2010, pp. 590, 592). —H. Graham et al. (2006, p. 6). —I. Chase et al. (2006, four-gene MP 
tree, pp. 69-70), ML tree in this study, Soltis et al. (2011, p. 7). 


ecological data were conducted under ML models (an 
independent model, I, and a dependent model, D) 
using the Markov chain Monte Carlo (MCMC) method 
under the Discrete option in BayesTraits 1.0. A 
likelihood score representing goodness of fit was 
calculated for the two models (Table 1, Part II). The 
likelihood ratio (LR) was determined by LR — 2 
(likelihood [dependent model]|-likelihood |indepen- 
dent model]) and was tested against a y? distribution 
with four degrees of freedom (following Pagel, 1994). 
Transition rates for correlated characters were 
estimated following Pagel and Meade (2006), with 
eight possible character transitions for each of the six 
pairs of traits (Fig. 2). 

A total of 5,000,000 generations with a burn-in 
period of 10,000 generations and a sampling 
frequency of 100 generations was used for both 
dependent and independent models, using a revers- 
ible-jump hyperprior model, exponential prior distri- 
bution, and rate variations ranging from 8 to 20. The 
difference in the harmonic means of the log- 
likelihoods between the two models was calculated 
twice. The log Bayes factor value (log BF) was 
calculated to indicate the degree of evolutionary 
correlation between two characters: with values of two 
lo five representing a positive correlation, greater 


than five a strong correlation, and greater than 10 a 
very strong correlation (Pagel & Meade, 2006). 

We then tested specific hypotheses of contingent 
evolution by restricting two of the rates to be equal 
following the method of Pagel (1994) and Schmitt et 
al. (2009) to test if the loss/gain of a palynological 
character is related to the absence/presence of an 
ecological character (Table 1, Part II). The two rate 
coefficients per trail were created from a model that 
allows the characters to evolve independently on the 
tree in the BayesDiscrete tests for correlated 
evolution. We compared this seven-parameter re- 
stricted dependent model to a full dependent model 
using the likelihood ratio test with one degree of 
freedom based on Pagel (1994). That is, a significant 
likelihood ratio would indicate that the parameters 
were significantly different from one another, dem- 
onstrating that the state of one trait influences the 
Finally, the 
exclusiveness of character state changes was tested 


direction of evolution of the other. 


by restricting one of the rates to zero and recalcu- 
lating the likelihood ratio of the model. We then 
compared the restricted seven-parameter model to the 
unrestricted dependent model using a x? test with one 
degree of freedom. A significant value would suggest 
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Rate parameters for the eight possible transitions between four pairs of binary pollen characters and binary 


ecological characters under a model of dependent evolution, exine (character [char.] al) and helo-hydrophytes (char. bl) or 
hydrophytes (char. b2), tectum (char. a2) and helo—hydrophytes (char. bl) or hydrophytes (char. b2), eutectum (char. a3) and 
helo-hydrophytes (char. bl) or hydrophytes (char. b2), infratectum (char. a4) and helo—hydrophytes (char. bl) or hydrophytes 
(char. b2); characters and states follow Table 1. Arrows indicate direction of transition (only one state transfer while another is 


stable). 


that a character trail was not gained or lost in the 
presence of another character trait. 

Three pollen binary characters, eutectum, colu- 
mellate infractectum, and granular infractectum, were 
not analyzed by the contingent evolution test because 
they were confirmed to be negatively correlated with 
helo-hydrophyte plants in the BayesDiscrete test. 


RESULTS 
PALYNOLOGICAL OBSERVATIONS 


The pollen morphology of 20 representative taxa 
observed under LM and SEM is described below to 
exemplify the diversity of pollen characters and their 
states in the monocots (Figs. 3-6). The species are 
presented in phylogenetic order, corresponding to 
their placement in Figure 7. 

Acorus gramineus Sol. ex Aiton (Acoraceae, 
Acorales, Fig. 3A-C) has pollen dispersed as 
monads; grains heteropolar, bilaterally symmetrical, 
boat-shaped, amb usually elliptic; size small (ca. 13— 
20 um); monosulcate, aperture distal, with sculptured 
aperture membrane; tectum with scattered perfora- 
lions. 

Epipremnum pinnatum (L.) Engl. (Araceae, Alis- 
matales, Fig. 3D—F) has pollen dispersed as monads; 
grains heteropolar, radially symmetrical, globose 
(spheroidal), amb circular; size medium (ca. 26-30 
um); zonate, aperture equatorial; tectum perforate to 
fossulate. 

Alisma orientale (Sam.) Juz. (Alismataceae, Alis- 
matales, Fig. 3G, H) has pollen dispersed as monads; 


grains apolar, asymmetrical, globose (spheroidal), 
amb circular; size small to medium (ca. 18-22 um); 
pantoporate, apertures global, with sculptured aper- 
ture membrane; supratectal elements echinate; 
tectum perforate. 

Zostera marina L. (Zosteraceae, Alismatales, Fig. 
3l, J) has pollen dispersed as monads; grains apolar, 
filiform; size gigantic (> 1000 um long); inaperturate; 
without exine. 

Petrosavia stellaris Becc. (Petrosaviaceae, Petrosa- 
viales, Fig. 3K, L) has pollen dispersed as monads; 
grains heteropolar, bilaterally symmetrical, boat- 
shaped, amb elliptic; size small (ca. 16-22 um); 
monosulcate, aperture distal; tectum reticulate. 

Dioscorea melanophyma Prain & Burkill (Dioscor- 
eaceae, Dioscoreales, Fig. 4A, B) has pollen 
dispersed as monads; grains subisopolar, bilaterally 
symmetrical, globose (oblate), amb elliptic; size small 
(ca. 14-17 pm); disulculate not situated at a pole, 
apertures equatorial; tectum fossulate to reticulate. 

Stemona tuberosa Lour. (Stemonaceae, Panda- 
nales, Fig. 4C, D) has pollen dispersed as monads; 
grains heteropolar, bilaterally symmetrical, boat- 
shaped, amb elliptic; size medium (ca. 26-30 um); 
monosulcate, aperture distal; tectum rugulate. 

Trillium maculatum Raf. (Melanthiaceae, Liliales, 
Fig. 4E, F) has pollen dispersed as monads; grains 
apolar, radially symmetrical, globose (spheroidal), 
amb circular; size medium (ca. 33-40 um); inapertu- 
rale; lectum rugulate. 

Tulipa iliensis Regel (Liliaceae, Liliales, Fig. 4G— 
I) has pollen dispersed as monads; grains heteropolar, 
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Figure 3. SEMs of monocot pollen, arranged following the topology of Chase et al. (2006). A—C. Acorus gramineus Sol. ex 
Aiton, Acoraceae, Acorales. —A. Equatorial view, showing monad, boat-shaped pollen grain, of small size, with monosulcate 
aperture. —B. Distal view. —C. Detail of tectum with scattered perforations. D-F. Epipremnum pinnatum (L.) Engl., Araceae, 
Alismatitales. —D. Equatorial view, showing expanded apertural ring (zonate pollen). —E. Polar view. —F. Detail of fossulate 
to reticulate tectum. G, H. Alisma orientale (Sam.) Juz., Alismataceae, Alismatales. —G. Whole grain, pantoporate pollen. —H. 
Detail of echinate supratectal elements. I, J. Zostera marina L., Zosteraceae, Alismatales. —I. Filiform, inaperturate pollen 
grain. —J. Detail of unsculptured, very thin, and highly reduced exine. K, L. Petrosavia stellaris Becc., Petrosaviaceae, 
Petrosaviales. —K. Polar view, showing boat-shaped pollen grain of small size. —L. Detail of reticulate tectum. 
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Figure 4. SEMs of monocot pollen, arranged following the topology of Chase et al. (2006). A, B. Dioscorea melanophyma 
Prain & Burkill, Dioscoraceae, Dioscorales. —A. Polar view, showing disulculate apertures. Note that the pollen grain shown 
here is rather anomalous; grains of most examples of this taxon are rather more globose. —B. Detail of fossulate to reticulate 
tectum. C, D. Stemona tuberosa Lour., Stemonaceae, Pandanales. —C. Equatorial view, showing monosulcate aperture. —D. 
Detail of rugulate tectum. E, F. Trillium maculatum Raf., Melanthiaceae, Liliales. —E. Spheroidal and inaperturate pollen 
grain. —F. Detail of rugulate tectum. G-L. Tulipa iliensis Regel, Liliaceae, Liliales. —G. Equatorial view, showing colpi located 
at distal and proximal poles. —H. Oblique equatorial view, showing two colpi located at the proximal pole. —1. Detail of 
rugulate tectum. J-L. Allium cyaneum Regel, Liliaceae, Liliales. —J. Equatorial view, pollen grain of medium size. —K. Distal 
polar view. —L. Detail of striate tectum. 
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Figure 5. SEMs of monocot pollen, arranged following the topology of Chase et al. (2006). A-C. Hemerocallis fulva (L.) L., 
Xanthorrhoeaceae, Asparagales. —A. Equatorial view, showing very large-size grain. —B. Distal polar view. —C. Detail of 
reticulate tectum. D—F. Calamus nambariensis Becc., Arecaceae, Arecales. —D. Equatorial view. —E. Polar view, showing two 
equatorially located colpi. —F. Detail of perforate to reticulate tectum. G, H. Calectasia cyanea R. Br., Dasypogonaceae, 
Dasypogonales. —G. Distal polar view, showing monosulcate, operculate aperture. —H. Detail of foveolate tectum. I, J. 
Anigozanthos viridis Endl., Haemodoraceae, Commelinales. —I. Polar view, showing two equatorially located pores. —J. Detail 
of rugulate tectum. K, L. Cartonema parviflorum Hassk., Commelinaceae, Commelinales. —K. Distal polar view, showing 
sculptured aperture membrane. —L. Detail of areolate to rugulate tectum. 
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Figure 6. SEMs of monocot pollen, arranged following the topology of Chase et al. (2006). A-C. Heliconia subulata Ruiz & 
Pav., Heliconiaceae, Zingiberales. —A. Distal polar view. —B. Heteropolar, oblate, large, inaperturate pollen grain. —C. 
Surface of distal polar view, showing granular exinous residues. D—F. Carex scaposa C. B. Clarke, Cyperaceae, Poales. —D. 
Equatorial view, showing single poroid aperture at the distal pole and four to five apertures in the equatorial region. —E. 
Proximal polar view. —F. Detail of gemmate supratectal elements. G, H. Juncus sikkimensis Hook. f., Juncaceae, Poales. —G. 
Entire tetrad. —H. Detail of areolate tectum and gemmate supratectal elements. I, J. Eriocaulon cristatum Mart., Eriocaulaceae, 
Poales. —I. Small, spiral pollen grain. —J. Detail of echinate supratectal elements. K, L. Dendrocalamopsis beecheyana (Munro) 
Keng f. var. pubescens (P. F. Li) Keng f., Poaceae, Poales. —K. Equatorial view, showing single aperture at the distal pole. —L. 
Detail of areolate tectum; arrows indicate pore with an operculum and annulus. 
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Figure 7. Inferences of character state change for 19 pollen characters and two ecological characters reconstructed with 
hierarchical Bayesian inference using comprehensive coding on a maximum likelihood phylogeny, based on the molecular data 
of Chase et al. (2006). Character state changes assume a single transition between differing states at adjacent nodes. Changes on 
terminal branches are not shown. Numbers on bars indicate the character state transition from the lower to upper node. 
Character boxes and the number within indicate character state change, following Table 1. $, bootstrap (BP) = 100; *, BP 90%- 
99%. —A. Seven orders (Acorales, Alismatales, Petrosaviales, Dioscoreales, Pandanales, Liliales, and Asparagales) and 


outgroups are shown. 
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Figure 7. Continued. —B. The commelinid clade is shown. 


bilaterally symmetrical, boat-shaped, amb elongate; 
size medium (ca. 40-45 um); trisulcate, one aperture 
distal, two apertures proximal, with sculptured 
aperture membrane; tectum rugulate. 

Allium cyaneum Regel (Amaryllidaceae, Aspar- 
agales, Fig. ALL) has pollen dispersed as monads; 
grains heteropolar, bilaterally symmetrical, boat- 
shaped, amb elliptic; size medium (ca. 28-40 um); 


monosulcate, aperture distal; lectum striate. 


remainder of tree 
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Hemerocallis fulva (L.) L. (Xanthorrhoeaceae. 


Asparagales, Fig. 5A—C) has pollen dispersed as 


monads; grains heteropolar, bilaterally symmetrical, 


boat-shaped, amb oblong; size large to very large (ca. 


85-105 um); monosulcate, aperture distal; tectum 


reticulate. 


Calamus nambariensis Becc. (Arecaceae, Are- 


cales, Fig. 5D—F) has pollen dispersed as monads; 


grains heteropolar, bilaterally symmetrical, globose 
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(spheroidal), amb elliptic; size medium (ca. 30-38 
um); disulcate, aperture equatorial; tectum perforate- 
reticulate. 

Calectasia cyanea R. Br. (Dasypogonaceae, Dasy- 
pogonales, Fig. 5G, H) has pollen dispersed as 
monads; grains heteropolar, bilaterally symmetrical, 
boat-shaped, amb oblong; size medium to large (ca. 
38-56 gm); monosulcate, aperture distal with a 
unique operculum, not reported elsewhere in angio- 
sperms, comprising a sulcus provided with three 
elongated, more or less parallel strips of membrane 
(Chanda & Ghosh, 1976); tectum foveolate-reticulate. 

Anigozanthos viridis Endl. (Haemodoraceae, Com- 
melinales, Fig. 51, J) has pollen dispersed as monads; 
grains subisopolar, bilaterally symmetrical, globose 
(oblate), amb elongate; size medium to large (ca. 45- 
65 gum); diporate, apertures equatorial; tectum 
rugulate. 

Cartonema parviflorum Hassk. (Commelinaceae, 
Commelinales, Fig. 5K, L) has pollen dispersed as 
monads; grains heteropolar, bilaterally symmetrical, 
boat-shaped, amb oblong; size medium (ca. 35-42 
um); monosulcate, aperture distal; tectum areolate to 
rugulate. 

Heliconia subulata Ruiz & Pav. (Heliconiaceae, 
Zingiberales, Fig. GA—C) has pollen dispersed as 
monads; grains heteropolar, radially symmetrical, 
globose (oblate), amb circular; size large (ca. 64—68 
um); inaperturate; exineless with granulate elements. 

Carex scaposa C. B. Clarke (Cyperaceae, Poales, 
Fig. 6D—F) has pollen grains dispersed as “pseudo- 
monads" (permanent tetrads with three rudimentary 
grains); grains heteropolar, asymmetrical, globose 
(prolate), amb circular; size medium (ca. 24-28 um); 
polyporate (usually with five or six poroids), apertures 
distal or equatorial, respectively; supratectal ele- 
ments gemmate; tectum perforate. 

Juncus sikkimensis Hook. f. (Juncaceae, Poales, 
Fig. 6G, H) has pollen dispersed as tetrads; tetrads 
radially symmetrical and globose (spheroidal); tetrads 
of medium size (ca. 34-38 um); individual grains 
heteropolar and radially symmetrical, globose (sphe- 
roidal), amb circular; monoporate, aperture distal; 
supratectal elements gemmate; tectum areolate with 
gemmae and perforations. 

Eriocaulon cristatum Mart. (Eriocaulaceae, Poales, 
Fig. 6L, J) has pollen dispersed as monads; grains 
apolar, asymmetrical, globose (spheroidal), amb 
circular; size small (ca. 16-18 um); aperture spiral; 
supratectal elements echinate; tectum perforate. 

Dendrocalamopsis beecheyana (Munro) Keng f. var. 
pubescens (P. F. Li) Keng f. (Poaceae, Poales, Fig. 6K, 
L) has pollen dispersed as monads; grains hetero- 
polar, radially symmetrical, globose (spheroidal), amb 
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circular; size medium (ca. 30—34 um); monoporate, 
operculate, with a surrounding annulus, aperture 
distal; supratectal elements gemmate; tectum areo- 
late. 


POLLEN VARIATION 


Most of the monocots studied possess heteropolar 
pollen grains dispersed as monads, with distal 
monosulcate or monosulcate-derived (e.g., a sulcus 
reduced to a pore) apertures, and a columellate 
infratectum, with an exine and foot layer but without 
an operculum or annulus. In addition to these 
generally consistent character states, apertural char- 
aclers in particular show taxonomically ordered 
patterns of variation: inaperturate pollen predomi- 
nates in the orders Alismatales and Zingiberales, 
whereas monosulcate pollen occurs in most taxa 
within the Liliales and Asparagales, in near-basal 
lineages such as Butomus L. (Butomaceae, Alisma- 
tales), in Gymnostachys R. Br. and Orontium L. 
(Araceae, Alismatales), and in some genera of 
Bromeliaceae (Poales). Highly specialized and 
unique pollen shapes (such as filiform in Zostera 
L.) with a highly reduced exine occur in the seagrass 
families (Cymodoceaceae, Potamogetonaceae, and 
Zosteraceae). Pollen morphological variation was 
observed to exist even within genera (in our 
comprehensive dataset for monocots, Appendix 2), 
with especially high variation found in Dioscorea L. 
(Dioscoreales), Pandanus Parkinson, and Stemona 
Lour. (both in Pandanales), as well as Trillium L. 
(Liliales), with lower variability seen in Alisma L. 
(Alismatales), Canna (Zingiberales), Juncus L., and 


Typha L. (both Poales). 


ANCESTRAL STATE RECONSTRUCTION 


The topology of the new ML tree based on DNA 
data from Chase et al. (2006) used for character 
optimization showed some incongruence with the 
original seven-gene tree derived from maximum 
parsimony (MP) analysis (Chase et al., 2006). 
Discrepancies included the positions of Boryaceae. 
Campynemataceae, Dasypogonaceae, Doryanthaceae, 
Heliconiaceae, Ixioliriaceae, Mayacaceae, Musaceae, 
and Typhaceae, for the clade comprising Eriocaula- 
ceae and Xyridaceae, for the position of the genera 
Aphelia R. Br., Baloskion Raf., and Elegia L., as well 
as internal relationships within the Asparagaceae s.l., 
which includes Agavaceae, Aphyllanthaceae, Hya- 
cinthaceae, Laxmanniaceae, Ruscaceae, and Themi- 
daceae (Angiosperm Phylogeny Group III, 2009). 
Nevertheless, these conflicts were not well supported, 
with mainly low (bootstrap values below 5096) or 
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moderate (bootstrap values 50%-80%) support in one 
or both trees (cf. Fig. 7). 


Of the 19 pollen characters 
studied, only pollen dispersal unit (cf. Table 1, 


Pollen characters. 


character |char.] 1) evolved unequivocally throughout 
all nodes with all five methods. Almost all pollen 
characters displayed unambiguous plesiomorphic 
states for monocots as a whole in all five methods 
(17 characters in analysis of the comprehensive 
dataset with Fitch parsimony inference [CFP], 19 in 
analysis of the comprehensive dataset with hierar- 
chical Bayesian inference [CHB], 17 in analysis of 
the democratic dataset with Fitch parsimony infer- 
ence [DFP], 19 in analysis of the democratic dataset 
with maximum likelihood inference [DML], and 18 in 
analysis of the democratic dataset with hierarchical 
Bayesian inference [DHB]). The most consistent 
plesiomorphic states across all five methods were 
dispersal unit (char. 1) as monad, polarity (char. 2) as 
heteropolar, pollen size (char. 7) as medium, aperture 
number (char. 8) as one, aperture membrane (char. 
10) as sculptured, aperture shape (char. 11) as 
colpate, annulus (char. 12) as absent, exine extent 
(char. 13) as fully present, supratectal elements (char. 
15) as absent, infratectum structure (char. 18) as 
columellate, and foot layer (char. 19) as present. 
Unambiguous (or probable) plesiomorphic states for 
the remaining seven characters were not consistently 
reconstructed across all five methods (char. 5, shape 
class, and char. 16, supratectal element shape, are 
inapplicable states for this node). For symmetry in 
polar view (char. 3), pollen grains were bilateral using 
CFP, CHB, DML, and DHB, and bilateral or radial 
under DFP. For basic pollen shape (char. 4), the 
plesiomorphic state was found to be boat shaped 
using CHB, DFP, DML, and DHB, and boat shaped 
or globose under CFP. For outline in polar view (char. 
6), pollen grains were found to be circular using DFP 
and DML, elliptic under CHB, oblong under DHB, 
and circular or elliptic under CFP. For aperture 
position (char. 9), a distal position was recovered 
using DFP, DML, CFP, and CHB, and an equatorial 
one under DHB; for tectum (char. 14), pollen was 
semitectate using CFP, DFP, DML, and DHB, and 
eutectate under CHB. For tectum sculpture (char. 
17), pollen was recovered as reticulate using CHB, 
DFP, DML, and DHB, and perforate or reticulate 
under CFP. 

All 19 pollen characters studied were inferred to 
exhibit unambiguous state changes on at least one 
lineage of monocots, using all five analytical 
methods. The state changes observed at family level 
or above are summarized in Table 2. Tectum 
sculpture, aperture membrane, and pollen size for 
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characters 17, 10, and 7, respectively, had the 
greatest number of inferred state changes. Pollen 
dispersal unit (char. 1), annulus (char. 12), infratec- 
tum structure (char. 18), and foot layer (char. 19) had 
the fewest number of changes. Of the 12 orders of 
monocols investigated, Zingiberales were subtended 
by the greatest. number of state changes (9), 
Alismatales and Asparagales by the fewest (zero). 


Ecological characters. The two ecological char- 
acters, pollination type and habitat moisture, were 
optimized onto the new ML tree using CFP and CHB 
methods. For pollination type (char. 20), the 
plesiomorphic state for monocots was found to be 
wind pollination (under both CFP and CHB). Water 
pollination was inferred to have evolved on the 
lineage leading to alismatids (CHB). Herein, alisma- 
tids we sampled include the Alismataceae, Butoma- 
ceae, Hydrocharitaceae, Cymodoceaceae, Juncagina- 
ceae, Potamogetonaceae, and Zosteraceae (Fig. 7). 
Shifts to wind pollination occurred in the clade 
comprising the alismatid families Juncaginaceae, 
Potamogetonaceae, and Zosteraceae (CHB), in the 
eraminid—restiid clade (CFP), in the Poales excluding 
Bromeliaceae and Typhaceae (CHB), Typhaceae 
(CFP, CHB), and Thurniaceae (CFP). Beyond the 
Poaceae, graminids include the Ecdeiocoleaceae, 
Joinvilleaceae, and Flagellariaceae; restiids encom- 
pass the Anarthriaceae, Centrolepidaceae, and Res- 
tionaceae. 

For habitat moisture (char. 21), the plesiomorphic 
state for monocots was inferred as mesic (CFP and 
CHB). Xeric or drier habitats evolved on the clade 
comprising Haemodoraceae and Pontederiaceae 
(CHB). Wetter habitats were found to evolve on six 
lineages including the clade comprising Commeli- 
nales and Zingiberales, and Thurniaceae (Poales; 
CFP or CHB). Wet habitats evolved on four lineages, 
i.e., Acorales, core alismatids (Alismatales), Typha- 
ceae (Poales), and the clade comprising Eriocaula- 


ceae and Xyridaceae (CFP and/or CHB; cf. Table 2). 


POLLEN CHARACTER STATE CHANGES 


The minimum number of state change steps of 
each of the 19 pollen characters optimizing on the 
nine alternative topologies tested (Fig. 1) were 
counted in the CFP analysis. Fifteen characters 
(excluding symmetry, basic shape, shape class, and 
outline in polar view [amb]) possessed the same 
minimum number of steps throughout all nine 
topologies. Within the topologies, A through F focus 
on the controversial relationships among major clades 
within commelinids, and topologies A, D, and E yield 
the fewest steps for symmetry and basic shape. 
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Table 2. Pollen character state changes inferred from comprehensive and democratic datasets of 19 pollen morphological and 


two ecological characters for lineages at or above family level in monocots analyzed with Fitch parsimony (FP), maximum likelihood 


(ML), and hierarchical Bayesian inference (HB). Pollen characters and their states (in brackets) correspond to Table 1. Superscripts 


refer to the number of the method used: 'CFP, ?CHB, *DFP, *DML, and °DHB. 


Taxon 


Acorales 
Alismataceae 


Alismatales 
Alismatids 
Alstroemeriaceae 
Amaryllidaceae 
Araceae 
Arecales 


Asparagaceae 

Asparagales 

Asparagales-commelinids 

Asparagales with Orchidaceae 
excluded 

Astelioids 

Boryaceae 

Bromeliaceae 

Colchicaceae 

Commelinaceae 

Commelinaceae-Hanguanaceae 

Commelinales 

Commelinales-Zingiberales 

Commelinales—Zingiberales—Poales 

Commelinids 

Commelinids with Arecales excluded 

Core asparagoid 

Cyclanthaceae 

Cyclanthaceae—Pandanaceae 


Cyperaceae 
Cyperids 


Dasypogonales 

Dioscoreaceae 

Dioscoreales 
Dioscoreales—Pandanales 
Eriocaulaceae—X yridaceae 
Graminids 
Haemodoraceae—Pontederiaceae 
Ixioliriaceae-Tecophilaeaceae 
Juncaceae 

Liliaceae 

Liliales 
Lowiaceae-Strelitziaceae 
Melanthiaceae 

Monocots 
Monocots—magnoliids 
Musaceae—Heliconiaceae 


Nartheciaceae 
Orchidaceae 
Pandanaceae 
Pandanales 


Characters and states 
Pollen characters 


7 (029, 14 (21%), 17 (1 2285) 


3 (23), 4 (13), 5 (13), 6 (03), 7 (03), 8 (428), 11 (1), 


14 (225), 15 (122), 17 (12) 

None 

2 (0255), 3 (12), 4 (19), 6 (03, 23), 7 (43), 8 (0255) 

7 (28), 17 (33) 

6 (23), 7 (2°) 

14 (2125), 17 (12255) 

2 (2255, 4 (14), 6 (22), 9 (225, 11 (22), 14 (23), 15 
13:17 (13) 


None 


None 

14 (13) 

6 (12), 7 (23) 
15 (135, 17 (13) 
7 (233) 


2 (02, 3 (125, 7 (05 , 13 (03) 

6 (13) 

7 (23, 33, 8 (02), 13 (03) 

S (125). 4 (155 5025 

4 (13), 10 (0%) 

6 (23), 7 (23), 9 (03), 14 (21) 

None 

None 

SE, 105.10 (05, TE 152), 14:275, 
17 (135 

3 (23) 


1 (13855, 6 (02), 7 (03), 11 (13), 14 
(03), 17 (02, 14) 

7 (24), 14 (2%), 17 (12) 

9 (22), 14 (238), 17 (13) 

5 (12), 9 23), 10 (13, 22, 16 (13) 

None 

2 (03), 6 (13), 14 (13) 

None 

2 (223, 7 (23, 32), 8 (2%, 32), 14 (23) 

6 (13, 10 (25, 14 (23), 17 (335) 

None 

6 (23) 

6 (23) 

Qa, 34) 

(025, 13 (03), 19 (03) 

(135), 7 (03), 9 (02), 10 (03), 14 (12) 

(23) 

None 

7 (0%), 10 (235, 17 (13) 

10 (233, 11 (13) 

5 (12, 23, 6 (23), 7 (028), 12 (13), 13 

14 (23) 


(24), 15 (14), 16 


7 
2 
6 
6 


(13), 17 (03) 


Ecological characters 


21 (33) 


20 (32), 21 (32) 


20 (53) 


20 (51) 


20 (42) 
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Table 2. Continued. 
Characters and states 
Taxon Pollen characters Ecological characters 
Petermanniaceae—A Istroemeriaceae— 7 (533) 
Colchicaceae 
Petrosaviales 7 (02%), 14 (03) 21 (23) 
Philesiaceae-Rhipogonaceae 2 (0), 3 (13), 4 (12), 6 (13), 8 (03), 15 (13, 17 (03) 
Philydraceae-Haemodoraceae— 8 (225), 13 (13 
Pontederiaceae 
Poaceae 10 (24), 16 (03), 17 (05) 
Poales 3 (15, 5 (13), 6 (05, 11 (13) 
Poales with Bromeliaceae and 11 (15, 14 (22), 15 (125), 17 (13) 2003) 
Typhaceae excluded 
Restiids 10 (25), 17 (13) 
Restiids—graminids 6 (03, 11 (18), 12 (14), 15 (14), 17 (14) 20(15) 


14 (03, 12) 


Thurniaceae ri 


Stemonaceae 
Tofieldiaceae 24), 8 (24) 
Typhaceae 
Velloziaceae 
g, 
Xyridaceae 
Zingiberaceae—Costaceae 


(0 
2 ( 
3 ( 
7 ( 
Xanthorrhoeaceae 6 (1 
3( 
3 ( 
2 ( 


Zingiberales 


4). 16 (13), 17 (02) 


13), 4 (13), 6 (03), 11 (1 ze 
025), 10 (22), 14 (12), 17 (33) 


20(15, 21(27) 


20(12), 21(32) 


2, 5 (22), 7 (22), 8 (023), 13 (0%), 19 (12) 
225). 7 (33), 8 (43), 13 (13) 
028), 3 (15, 4 (13, 


23), 6 (03), 7 (23, 8 (023), 13 


(0 134) 14 (2 3), 19 (0 12345) 


Topology G possesses the fewest steps for basic 
shape, shape class, and outline in polar view (amb), 
of the three topologies (G—I) that focus on the 
basalmost branching group of Poales. 

To illustrate character state changes throughout the 
phylogeny of monocots, we mapped the unambiguous 
character state changes in the 19 pollen characters 
and two ecological characters (Table 1) at internal 
nodes based on CHB optimization (Fig. 7), the 
method that tends to produce the fewest equivocal 
states at internal nodes in our analysis (Fig. 7). 
Changes in pollen character states were not distrib- 
uted regularly but were clustered in certain lineages. 
While most nodes were associated with five or fewer 
changes, 14 internal nodes displayed at least five 
inferred changes; 10 nodes were seen with five 
changes, one node with six changes, two nodes with 
seven changes, and one node with nine changes. The 
10 nodes with five changes include the root nodes of 
alismatids, the Stemonaceae, Commelinales, Zingi- 
berales, and Typhaceae. The one node with six 
changes was the root of Xyridaceae. The two nodes 
with seven changes were the clade comprising 
Eriocaulaceae and Xyridaceae, and the Arecales. 
The root node of Alismataceae displayed nine 
changes. 

Five orders displayed more than 20 internal pollen 
character state changes: Alismatales (41 changes), 


Asparagales (24), Liliales (21), Pandanales (22), and 


Poales (45). In contrast, the Dioscoreales, Arecales, 
and Dasypogonales possessed very few or no changes: 
four, three, and zero, respectively. Taking the relative 
sampling and taxonomic size of a lineage into account 
(pollen character state changes per internal node 
within a clade), the highest relative rate of pollen 
character state changes was found in Alismatales, 
and Arecales (3.73, 3.25, and 3, 


respectively), while Zingiberales, Asparagales, and 


Commelinales, 


Dasypogonales had the lowest relative rate of pollen 
character change (1.17, 0.83, and 0, respectively). 

Thus, these 19 pollen characters (Table 1) showed 
variable frequencies of state change during the 
evolution of monocots. Changes of character state at 
internal nodes were found to be more frequent in 
pollen size (char. 7, with 29 changes), pollen outline 
in polar view (char. 6, 22 changes), and tectum 
sculpture (char. 17, 23 changes), whereas pollen 
dispersal unit (char. 1, one change), annulus (char. 
12, two changes), infratectum structure (char. 18, one 
change), and pollen foot layer (char. 19, two changes) 
had far fewer (Fig. 7). 


ANALYSES OF CORRELATED EVOLUTION 


With the exception of exine presence (MCMC 
method, log Bayes factor value [log BF] = 3.19), we 
failed to find positive correlations between exine 
stratification states and hydrophytes using ML and 
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Results of correlated evolution analyses using maximum likelihood (ML) and Markov chain Monte Carlo (MCMC) 


methods. Characters and character states are numbered, following Table 1. 


Independent 


Trait pair (pollen and ecological character) model 


Exine stratification and helo-hydrophytes 
Exine (character [char.] al) and helo- 


hydrophytes (char. bl) —94.01 
Tectum (char. a2) and helo—hydrophytes 

(char. bl) —102.50 
Eutectum (char. a3) and helo—hydrophytes 

(char. bl) —146.43 
Infratectum (char. a4) and helo—hydrophytes 

(char. bl) —98.90 
Columellae (char. a5) and helo—hydrophytes 

(char. bl) —109.91 
Granules (char. a6) and helo-hydrophytes 

(char. bl) —82.80 

Exine stratification and hydrophytes 

Exine (char. b1) and hydrophytes (char. b2) —68.09 
Tectum (char. a2) and hydrophytes (char. 

b2) —76.58 
Eutectum (char. a3) and hydrophytes (char. 

b2) —120.51 
Infratectum (char. a4) and hydrophytes 

(char. b2) —72.97 
Columellae (char. a5) and hydrophytes 

(char. b2) —83.99 
Granules (char. a6) and hydrophytes (char. 

b2) —56.88 


* P < 0.05, LBF > 2. 
* P < 0.01, LBF > 5. 
t P < 0.005, LBF > 10. 


MCMC analysis methods (Table 3). ML and MCMC 
analyses of correlation for the eutectum and helo— 
hydrophily yielded inconsistent results (P > 0.05, log 
BF = 2.85, Table 3). 

There was very strong or strong support for 
correlated evolution between exine presence and 
habitat for helo-hydrophytes (both with significance 
of P < 0.01 and log BF > 10), and strong or positive 
support for correlated evolution between three exine 
characters (presence of tectum, eutectum, and 
infratectum) and helo—hydrophily (P < 0.01, log 
BF = 5-10). The remaining two pollen traits tested 
(columellate infratectum and granular infractectum) 
were nol significantly correlated with helo-hydro- 
phily (a negative correlation with P > 0.05 and log 
BF < 2). Results of correlated evolution analyses for 
the seven pairs of traits are shown in Table 3. 

In a test of specific hypotheses with regard to 
contingent evolution, significant support was found: 
the presence of pollen exine was not gained 
preferentially in helo-hydrophytes compared to those 


Likelihood 
Dependent likelihood 


ML MCMC 
Harmonic mean of log-likelihood 


Independent Dependent Log Bayes 


model ratio (LR) model model factor (LBF) 
—87.12 13.77 99:21 —92.71 13.12? 
—95.33 14.33! —106.99 — 102.53 8.931 
—142.33 8.22 —150.63 —149.21 2.85* 
—92.20 15.40! —104.04 —101.73 4.62* 
—105.22 9.38 113.62 113.72 —0.20 
82.33 0.93 —89.06 —88.87 0.39 
—63.61 8.96 =(1,20 —69.60 3.19* 
—193.15 6.86 — 718.49 — 718.20 0.57 
—118.61 3.02 — 126.46 —125.69 1.54 
—70.11 9.14 —(5.15 —15.33 0.83 
—81.69 4.60 —87.19 —87.59 —0.81 
—56.63 0.49 —61.97 —62.19 —0.65 


that were not helo—hydrophytes (hypothesis q24 < 
q13, probability values 0.02). Meanwhile, in a test of 
specific hypotheses of exclusiveness evolution, the 
loss of pollen exine was not significant in plants that 
were not helo-hydrophytes (hypothesis q31 — 0, 
probability values 0.1304) but was significant in 
helo—hydrophytes (hypothesis q42 — 0, probability 
values 0.0120). The losses of pollen tectum and 
infratectum were significant in both helo—hydro- 
phytes and those that were not helo—hydrophytes 
(hypothesis q31 = 0 rejected, probability values 
0.0347 and 0.0262 respectively, P < 0.005). 


DISCUSSION 


SYSTEMATIC SIGNIFICANCE OF POLLEN MORPHOLOGY IN 
MONOCOTS 


Our optimizations using five methods (analysis of a 
comprehensive dataset with Fitch parsimony and 
hierarchical Bayesian inference [abbreviated as CFP 
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and CHB] and analysis of a democratic dataset with 
Fitch parsimony, maximum likelihood, and hierar- 
chical Bayesian inference [abbreviated as DFP, 
DML, and DHB respectively]) allowed us to infer 
the pollen morphological state changes or synapo- 
morphies for multiple branches within the monocot 
phylogenetic tree. We found that pollen characters 
may be of systematic value for delimiting lineages 
and may provide clues to controversial relationships 
within monocots (cf. Table 2). Each of the major 
lineages of monocots is discussed below. 


Acorales. The order Acorales was recognized in 
both Angiosperm Phylogeny Group II (2003) and 
Angiosperm Phylogeny Group III (2009) classifica- 
tions. Acorus, the single genus in the order, was 
previously placed in Araceae (reviewed in Grayum, 
1987) but later found to be sister to the rest of the 
monocots based on molecular data (e.g., Chase et al., 
1993; Duvall et al., 1993a, 1993b; Chase et al., 
2006; Graham et al., 2006). Its isolated position is 
strongly congruent with several synapomorphies, 
including small size (char. 7, CFP, CHB, and 
DML), eutectate pollen (char. 14, CFP, DFP, DML, 
and DHB), and perforate tectum (char. 17, CHB, 
DFP, DML, and DHB) (Table 2). 


Alismatales. The large order Alismatales com- 
prises ca. 4560 herbaceous species of aquatic and 
moist habitats and is the second-diverging clade 
within the monocots (Angiosperm Phylogeny Group 
III, 2009). No unambiguous pollen synapomorphies 
were found for the Alismatales. Recent molecular 
phylogenetic estimates have identified three major 
clades within Alismatales: Araceae, Tofieldiaceae, 
and the alismatid aquatic. families (Chen et al., 
2004a, 2004b; Chase et al., 2006), all of which are 
characterized by character state changes. For Ara- 
ceae, states include pollen eutectate (char. 14, CFP, 
DFP, DML, and DHB) with perforate tectum (char. 
17, CHB, DFP, DML, and DHB); Tofieldiaceae are 
characterized by isopolar and diaperturate pollen 
(char. 2 and char. 8, DML). Character states 
identified for the aquatic alismatid families are 
pollen apolar (char. 2, CHB, DML, and DHB), 
radially symmetrical in polar view (char. 3, CHB), 
elobose (char. 4, CHB and DHB), circular in polar 
view (char. 6, DHB), gigantic in size (char. 7, CHB 
and DHB), and lacking apertures (char. 8, CHB, 
DML, and DHB; Table 2). 


Petrosaviales. The family Petrosaviaceae (Kubitz- 
ki, 1998a) was newly recognized as the order 
Petrosaviales in Angiosperm Phylogeny Group IM 


(2009). Most members of the group were previously 
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treated as belonging to a diversity of families 
including Liliaceae (Ohba, 1984), Nartheciaceae, 
Tofieldiaceae, and Melianthaceae (Dahlgren et al., 
1985). The isolated position of this order is 
characterized by two pollen character state changes, 
to small pollen size (char. 7, CHB, DML, and DHB) 
and intectate pollen (char. 14, DHB; Table 2). 


Dioscoreales and Pandanales. Dioscoreales are 
moderately supported as sister to Pandanales based 
on molecular evidence (Chase et al., 2006; Graham 
et al., 2006; Petersen et al., 2006). We found no 
palynological characters to support this relationship. 
The Dioscoreales themselves were characterized by 
three pollen character state changes, to spheroidal 
pollen (char. 5, DHB), equatorial apertures (char. 9, 
CHB and DHB), and echinate supratectal elements 
(char. 16, CHB). The monophyly of the Dioscor- 
eaceae s.l. (including the former Taccaceae and 
Trichopodaceae, as accepted by the Angiosperm 
Phylogeny Group H based on the combined 
molecular-morphological analyses of Caddick et al. 
[2000, 2002a, 2002b]) is congruent with state 
changes to equatorial apertures (char. 9, CHB), 
eulectate pollen (char. 14, DFP, DML, and DHB), 
and perforate tectum sculpture (char. 17, DHB). 
Pollen morphology is particularly congruent with 
known systematic relationships within Pandanales. 
The sister relationship between Cyclanthaceae and 
Pandanaceae is congruent with pollen character state 
transitions to radial symmetry (char. 3, CFP and 
DML), globose shape (char. 4, CFP, CHB, DFP, and 
DML), oblate pollen (char. 5, CHB), smooth aperture 
membrane (char. 10, CFP, DFP, and DML), porate 
apertures (char. 11, CFP, DFP, and DML), eutectate 
pollen (char. 14, DFP and DML), and a perforate 
tectum (char. 14, DFP and DML). Pandanaceae as a 
group are defined by oblong outline in polar view 
(char. 6, CHB), small size (char. 7, DFP, DML, and 
DHB), annulus present (char. 12, CHB and DHB), 
atectate pollen (char. 13, DHB), and imperforate 
tectum sculpture (char. 17, DHB). Stemonaceae 
received no defining pollen character states, while 
Velloziaceae are characterized by small pollen size 
(char. 7, all five methods) with an operculum (char. 
10, CHB and DHB), semitectate (char. 14, CHB), 
and rugulate tectum (char. 17, CHB; Table 2). 


Liliales. 
to oblong outline in polar view (char. 6, CFP and 
DFP), is found on the stem of Liliales. The recent 
merger of Luzuriagaceae into Alstroemeriaceae 


Only one pollen character state change, 


(Mabberley, 2008; Angiosperm Phylogeny Group 
III, 2009) is congruent with both large pollen size 
(char. 7, DHB) and rugulate tectum (char. 17, DHB). 
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The placement of Uvulariaceae within Colchicaceae 
(Fay et al., 2006) is congruent with supratectal 
elements present (char. 15, CHB and DHB) and 
perforate tectum sculpture (char. 17, DHB). The 
placement of the Petermanniaceae has been debated 
for some time (Angiosperm Phylogeny Group ll, 
2003; Janssen & Bremer, 2004; Fay et al., 2006). 
The topology we used in our study, in which it is 
sister to the clade comprising Alstroemeriaceae and 
Colchicaceae (Fay et al., 2006), is congruent with two 
pollen character state transitions, to large size (char. 
7, DHB) and fossulate tectum (char. 17, CHB and 
DHB; Table 2). 


Asparagales. The placement of Asparagales as 


sister to commelinids has been moderately supported 
by Chase et al. (2006) and strongly by Graham et al. 
(2006). In our study, this relationship is congruent 
with only one pollen state (char. 10, smooth aperture 
membrane, DHB). No state change was found to 
occur on the branch leading to Asparagales. The 
monophyly of the astelioid clade (including Astelia- 
ceae, Blandfordiaceae, Boryaceae, Hypoxidaceae, 
and Lanariaceae) has been well supported in recent 
phylogenetic studies (Seberg et al., 2012; Chen el al., 
2013); however, no pollen character state changes 
were found for this clade. The placement of 
Ixioliriaceae as sister to Tecophilaeaceae was strongly 
supported in the study of Chen et al. (2013) and is 
congruent with transitions to elliptic pollen outline in 
polar view (char. 6, DHB), operculate aperture (char. 
10, DML), eutectate pollen (char. 14, DHB), and 
rugulate tectum (char. 17, CHB and DHB; Table 2). 


Commelinids. The commelinid clade is strongly 


supported in molecular studies (reviewed in Chase et 
al., 2006), and the monophyly of each of the five 
commelinid orders is well established (e.g., Davis et 
al., 2004; Chase et al., 2006). However, internal 
relationships among its five major groups remain 
controversial, particularly the positions of Dasypogo- 
nales and Arecales (e.g., Chase et al., 2006; Graham 
et al., 2006; Givnish et al., 2006, 2010; Angiosperm 
Phylogeny Group III, 2009; Soltis et al., 2011; Barrett 
et al., 2012; Davis et al., 2012). Six main conflicting 
topologies within the commelinids are summarized in 
Figure 1. 

Among the six alternative topologies tested (A—F), 
topologies A, D, and E produced the minimum 
number of pollen character state change steps in the 
CFP analysis, and, thus, they are inferred to be 
favored by pollen morphology based on the parsimony 
theory, i.e., (Arecales (Dasypogonales (Commeli- 
nales—Zingiberales, Poales))), ((Arecales, Dasypogo- 
nales) (Commelinales-Zingiberales, Poales)), and 
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((Arecales, Poales) (Commelinales—Zingiberales, Da- 
sypogonales)), respectively (Fig. 1). It is noteworthy 
that the latest phylogenetic studies using plastid 
genes and genomes (Givnish et al., 2010; Barrett et 
al., 2012; Davis et al., 2012) support Arecales plus 
Dasypogonales as sister to the rest of the commelinid 
monocots (topology D). 

The sister relationship between Commelinales and 
Zingiberales has been well supported (Chase et al., 
2006; Givnish et al., 2006, 2010; Graham et al., 
2006; Soltis et al., 2007, 2011; Qiu et al., 2010) and 
is also congruent with transitions to very large or large 
size (char. 7, CHB and DHB, respectively), inapertu- 
rate (char. 8, CHB) and exineless pollen (char. 13, 
CHB and DHB; Table 2). Both topologies A and D 
contain a clade comprising Commelinales and 
Zingiberales as sister to Poales, which is character- 
ized by radial symmetry (char. 3) in both. 

The studies of both Givnish et al. (2006) and 
Chase et al. (2006) presented consistent sister 
relationships between Haemodoraceae and Ponteder- 
iaceae. The clade comprising Haemodoraceae and 
Pontederiaceae is congruent with multiple pollen 
character state changes, to isopolar pollen (char. 2, 
CHB and DHB), very large size or large size (char. 7, 
CHB and DHB, respectively), two apertures (char. 8, 
DFP and DML), and eutectate grains (char. 13, DHB; 
Tables 1, 2). 

The Zingiberales are a well-characterized clade 
with nine inferred shared pollen character states: 
apolar (char. 2, CHB, DML, and DHB), radially 
symmetrical (char. 3, DML), globose pollen (char. 4, 
CHB), circular outline in polar view (char. 6, CHB 
and DHB), large size (char. 7, DML), no apertures 
(char. 8, CFP, DFP, DML, and DHB), no exine (char. 
13, CFP, DFP, and DML), eutectate pollen (char. 14, 
DHB), and no foot layer (char. 19, all five methods). 
Two pairs of sister families are found within the order: 
Costaceae and Zingiberaceae, and Lowiaceae and 
Strelitziaceae (Kress et al., 2001; Chase et al., 2006; 
Givinish et al., 2006). The former pairing is 
congruent with transitions to asymmetrical pollen 
(char. 3, CHB and DHB), very large size (char. 7, 
DHB), more than three apertures (char. 8, CFP, DFP, 
DML, and DHB), exineless (char. 13, CFP, DFP, and 
DML) and eutectate pollen (char. 14, DHB), and 
absence of the foot layer (char. 19, all five methods); 
the latter pairing is congruent only with very large 
pollen size (char. 7, DML and DHB). The monophy- 
letic group including Cannaceae, Marantaceae, 
Costaceae, and Zingiberaceae among other Zingiber- 
ales is congruent with one pollen character state 
change, to large pollen size (char. 7, CHB; cf. Fig. 1, 
Tables 1, 2). 
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Within Poales, the basalmost branching group has 
variously been recovered as Bromeliaceae, Typha- 
ceae, or Bromeliaceae-Typhaceae (Davis et al., 2004; 
Chase et al., 2006; Givnish et al 2006, 2010; 
Graham et al., 2006). These three conflicting basal 
topologies within Poales are shown in Figure 1G-I. 
Of these, topology G, with Bromeliaceae as the first 
lineage to diverge, and Typhaceae sister to the 
remaining Poales (including Rapateaceae and the 
eraminid, restiid, cyperid, and xyrid clades) requires 
the fewest pollen character state changes (cf. Fig. 
1G-L Table 1). This topology, which is consistent 
with the phylogenetic work of Givnish et al. (2010), is 
inferred to be favored by pollen morphology. 

The cyperid clade (including Cyperaceae, Junca- 
ceae, and Thurniaceae) has been well supported by 
molecular data (Chase et al., 2006; Givnish et al., 
2010) and can be characterized by pollen states 
including grains borne in tetrads (char. 1, all five 
methods), pollen circular outline in polar view (char. 
6, CHB), small grain size (char. 7, CHB and DHB), 
porate aperture (char. 11, CFP and DFP), supratectal 
elements present (char. 15, DML), gemmate supra- 
tectal element shape (char. 16, CHB and DHB), and 
imperforate tectum (char. 17, CHB and DHB). 
Herein, the pseudomonads of Cyperaceae are treated 
as hidden tetrads. The xyrids (including Eriocaula- 
ceae, Mayacaceae, and Xyridaceae) did not form a 
monophyletic group in Chase et al. (2006) and our 
study (Fig. 7). We found that a sister relationship 
between Eriocaulaceae and Xyridaceae would be 
characterized by apolar pollen (char. 2, CHB and 
DHD,, elliptic outline in polar view (char. 6, DHB), 
and semitectate pollen (char. 14, CHB). Mayacaceae 
can be distinguished from the other two families by 
its monosulcate and tectate to intectate pollen 
(Zavada, 1983), while there is no pollen character 
linking it with Eriocaulaceae and Xyridaceae. A 
sister relationship between the restiid and graminid 
clades is found to be congruent with pollen character 
state changes to circular outline in polar view (char. 
6, CHB), porate aperture (char. 11, CFP and DFP) 
with an annulus (char. 12, all five methods), 
supratectal elements present (char. 15, DML), and 
perforate tectum sculpture (char. 17, DML; cf. Fig. 7, 
Tables 1, 2). 


THE RELATIONSHIP BETWEEN MONOCOTS AND BASALLY 
BRANCHING ANGIOSPERMS 


The Chloranthaceae, magnoliids, monocots, Cera- 
tophyllaceae, and eudicots are all well supported as 
major lineages within angiosperms, but the relation- 
ships among them remain unclear (e.g., Qiu et al., 


2005, 2006, 2010; Soltis et al., 2005, 2008, 2011; 
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Moore et al., 2007). Studies have suggested several 
alternatives for the closest relatives of monocots, 
ranging from magnoliids (Hilu et al., 2003; Duvall et 
al., 2006) to the aquatic Ceratophyllaceae (Qiu et al., 
1999, 2005; Zanis et al., 2002), Ceratophyllaceae— 
eudicots (Moore et al., 2007; Soltis et al., 2011), 
Piperales (Duvall, 2000), or Laurales (Parkinson et 
al., 1999). In our study, although a sister relationship 
between monocots and magnoliids is congruent with 
shared pollen character states including oblong 
oulline in polar view (char. 6, DHB; cf. Fig. 7, Table 
1), an affinity of monocots to basal eudicots or 
Ceratophyllaceae has also been supported by pollen 
characters (Lu et al., 2015). 

Our studies suggest that convergent evolution in 
pollen character states occurs commonly between 
some monocols and basally branching angiosperms. 
Both Lu et al. (2015) and the present study suggest 
that both groups contain taxa with inaperturate pollen 
grains with a highly reduced exine and thick 
channeled intine (Zingiberales and Lauraceae; Stone, 
1987), trichotomosulcate pollen (Iridaceae, Areca- 
ceae, and Chloranthaceae; Chapman, 1987; Goldblatt 
& Le Thomas. 1992, 1997; Harley, 1999) or pollen 
wilh zonate aperlures (Araceae, Arecaceae, Nym- 
phaeaceae, and Eupomatiaceae; Woodland & Gar- 
lick, 1982; Grayum, 1992; Gabarayeva & Rowley, 
1994; Hesse et al, 2001; Harley € Dransfield, 
2003). Polyplicate striae are found in Araceae, 
Zingiberaceae, and Lauraceae (Hesse et al., 2000), 
and loss or reduction of exine is also common in both 
monocots and basal angiosperms. 


EVOLUTION OF POLLEN IN MONOCOTS 


In this study, plesiomorphic pollen morphological 
states for monocots were unambiguously inferred 
under all five optimization methods to include pollen 
dispersal as monads (char. 1, Table 1, Fig. 7), pollen 
heteropolarity (char. 2), medium size (char. 7), 
monosulcate apertures (characters 8 and 11) with a 
sculptured aperture membrane (char. 10) and lacking 
an annulus (char. 12), exine fully present (char. 13), 
supratectal elements absent (char. 15), infratectum 
structure columellate (char. 18), and foot layer 
present (char. 19). The earliest monocot pollen is 
therefore likely to have been characterized by most of 
these pollen features. 

Distal monosulcate pollen has long been consid- 
ered an ancestral state for both basal angiosperms 
and monocots (Wodehouse, 1936; Takhtajan, 1969; 
Zavada, 1983; Sampson, 2000; Doyle, 2005; Lu et 
al., 2015). This character is also shown unambigu- 
ously to be plesiomorphic in monocots under four of 
our five methods (the exception being DHB, which 
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inferred an equatorial aperture position to be 
plesiomorphic), and this is in agreement with 
previous studies (Penet et al., 2005; Ressayre el 
al., 2005). 

Evolution of the pollen dispersal unit (char. 1, 
Table 1, Fig. 7), infratectum structure (char. 18), and 
foot layer presence (char. 19) has been relatively 
conservative throughout the angiosperms (Worlley el 
al., 2015). These three pollen characters were found 
to be similarly stable in the evolution of monocots, 
having the lowest frequency of character state 
changes, despite the Orchidaceae (although here 
with only three taxa sampled) being noted as most 
diverse for pollen dispersal unit among angiosperms 
(Pacini & Hesse, 2002). Thus, the states dispersal as 
monads, columellate infratectum structure, and 
presence of a foot layer are found in the overwhelm- 
ing majority of monocots. Other character states are 
less commonly reported (e.g., granulate infratectum) 
or are restricted to certain derived, monophyletic 
groups (e.g., absence of the foot layer in Haemodor- 
aceae; Simpson, 1983). In contrast, pollen size (char. 
7) and tectum sculpture (char. 17) displayed greater 
frequency of change among monocots, as has also 
been the case throughout angiosperms (Wortley et al., 
2015). 

Pollen aperture characters have displayed frequent 
state changes throughout the phylogeny of monocots. 
Although the range of aperture types in monocots is 
fewer than that observed for other angiosperms 
(Harley & Zavada, 2000), apertures are key to the 
evolution of monocot pollen morphology. Two 
aperture character state shifts, from a proximal to a 
distal position in early seed plant evolution and to an 
equatorial position seen at the base of the eudicot 
clade, have been regarded as some of the most 
striking innovations in the evolutionary history of 
seed plants (Muller, 1970; Walker & Doyle, 1975; 
Chanda & Ghosh, 1979; Zavada, 1983, 1984a; 
Rudall & Bateman, 2007). There is no intermediate 
position of apertures between distal and equatorial, 
and the switch in position reflects developmental 
changes expressed during meiosis (Blackmore & 
Crane, 1998). In our analyses for monocots, under 
CHB, a further transition from distal to equatorial 
apertures is inferred at the root of the Arecaceae. 
Pollen grains ascribed to Arecaceae are abundant 
and widely distributed from the Late Cretaceous 
(Burnham & Johnson, 2004; Nichols & Johnson, 
2008) and include all four types of equatorial 
aperture found among monocots, from disulcate, 
found, for example, in Tofieldiaceae, Dioscoreaceae, 
and Liliaceae; to diporate, as seen in Bromeliaceae 
and Haemodoraceae; zonate, observed in Araceae, 
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Iridaceae, and Rapateaceae (Hesse & Zetter, 2005); 
and triaperturate, observed in Iridaceae, Commelina- 
ceae, and Orchidaceae (Harley, 2004). Further 
extensive sampling and analysis of pollen morphology 
of the Arecaceae and related taxa will provide 
valuable clues to our understanding of distal/ 
equatorial aperture transitions within angiosperms. 
Zavada (1983) identified a number of evolutionary 
trends in aperture type in the monocots, such as the 
loss of a sulcus and an increase in the number and 
types of apertures. The trend of loss of sulcus was 
further exemplified by both Furness and Rudall 
(2006) and our CHB analysis. Furness and Rudall 
(2006) revealed that there are apparently homoplastic 
trends in each family toward either reduction in 
aperture size (e.g., from sulci to pores) or to coverage 
of the aperture with exine (e.g., from monosulcate to 
monosulcate-operculate) within Pandanales. In our 
analysis, monosulcate apertures were inferred unam- 
biguously as plesiomorphic for monocots, whereas 
monoporate and inaperturate types were derived 
within the group. Monoporate pollen (char. 8 and 
char. 11) is inferred to have evolved on the branches 
leading to Cyclanthaceae—Pandanaceae (CFP, DFP, 
and DML), Typhaceae (all five methods), cyperids 
(CFP and DFP), and the graminid-restiid clade (CFP 
and DFP; Table 2). Within the clade comprising 
Philydraceae, Haemodoraceae, and Pontederiaceae, 
monosulcate pollen reoccurs in all taxa of the basally 
placed family Philydraceae (Simpson, 1985a), with 
diaperturate (char. 8, DFP and DML) or triaperturate 
pollen (char. 8, CHB) evolving in the more derived 
Haemodoraceae—Pontederiaceae (Table 2). Alismata- 
ceae, among the later-diverging lineages of Alisma- 
tales, are characterized by pantoporate pollen (char. 
8, CFP, CHB, DFP, and DML). Inaperturate pollen 
(char. 8, Fig. 7, Table 1) evolved more frequently on 
lineages with helo—-hydrophytic plants, such as 
alismatids (CHB, DML, and DHB) and within 
Zingiberales (CFP, DFP, DML, and DHB). 


ECOLOGICAL ADAPTATIONS IN POLLEN MORPHOLOGY 


Monocots display a high diversity of pollination 
types, including wind pollination (Poaceae and 
Cyperaceae, exceptionally insect pollination in 
bamboo grasses; Soderstrom & Calderon, 1971), 
water pollination (diverse seagrass lineages in 
Alismatales), and bird pollination (widespread across 
the eight families of the Zingiberales; Nur, 1976), as 
well as dispersal of aggregated pollen, usually by 
animal vectors (Orchidaceae; Harder, 2000). Each of 
these pollination types corresponds to a series of 
distinct pollen features. For example, bird-pollinated 
Zingiberales possess exineless pollen with a thick- 
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ened, elaborated intine to preserve water (Kress & 
Stone, 1982; Theilade & Theilade, 1996); the water- 
pollinated seagrasses in Alismatales have filiform and 
exineless pollen (Ackerman, 1995, 2000); and the 
wind-pollinated graminid clade of Poales have 
typically distally monoporate pollen with an opercu- 
lum and annulus (Linder & Rudall, 2005). It is 
interesting to note that the entomophilous bamboo 
genus, Pariana Aubl. has a unique pollen type 
without an annulus, which may imply a further 
possible association between pollen morphology and 
pollination type (Skvarla et al., 2003). In our CHB 
analysis, four internal nodes having not less than five 
changes in pollen character state were found to be 
closely associated with pollinator shifts: the root of 
the alismatids (from entomophily to hydrophily), the 
root of the Arecales (a reversal from ornithophily to 
entomophily), the root of Typhaceae (from ornithoph- 
ily to anemophily), and the root of the clade 
comprising Eriocaulaceae and Xyridaceae (a reversal 
from anemophily to entomophily; see Fig. 7). The 
interesting results here indicate that the transitions in 
pollen morphological states may be potentially 
correlated with the shifts in pollination type. 
Asparagales have the second highest diversifica- 
lion rate in monocots (Magallón & Castillo, 2009) yet 
the lowest average frequency of pollen character state 
changes at internal nodes among all orders except 
Dasypogonales (0.76). We found only 22 state 
changes across all 19 pollen characters studied in 
this clade: seven in pollen size (char. 7), four in 
outline in polar view (char. 6), three in aperture 
membrane (char. 6), four in tectum (char. 14), three 
in tectum sculpture (char. 17), and one in symmetry 
(char. 3), indicating that Asparagales pollen mor- 
phology has been consistent and comparatively stable 
during evolution (Fig. 7). In terms of ecological 
characters, habitat moisture (char. 21) was found to 
change frequently (seven times) within the core 
asparagoids, whereas pollination type (char. 20) did 
not change (Fig. 7): entomophily is consistently the 
most probable state under CHB at all internal nodes 
within Asparagales excluding the astelioids and 
Orchidaceae (Fig. 7). These observations suggest 
that changes in ecological interactions, particularly 
pollination type, may have remarkable impact as 
drivers of pollen morphological change within 
monocots, while pollen characters can remain 
remarkably stable under changing habitat moisture. 
Fourteen internal nodes in the phylogeny of 
monocots were found to be characterized by a high 
number of pollen morphological state changes under 
CHB (Fig. 7, Table 1), including the root node of 


Alismataceae (nine changes) the root nodes of 
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Arecales (seven changes) the clade comprising 
Eriocaulaceae and Xyridaceae (seven changes); 
Xyridaceae (six changes); 10 further nodes have five 
changes, including the roots of alismatids, Stemona- 
ceae, Commelinales, Zingiberales, and Typhaceae. 
Based on divergence time data for monocots provided 
by recent studies (Janssen & Bremer, 2004; Merckx 
et al., 2008; Bell et al., 2010; Smith et al., 2010; 
Bouchenak-Khelladi et al., 2014; Magallón et al., 
2015), the crown or stem node ages of 10 of these 14 
major clades are all dated to between 90 and 114 
million years ago, approximately corresponding to the 
Albian-Turonian stage of the Middle Cretaceous. 
Although many monocol orders are unrepresented by 
fossil records during the Late Cretaceous, there is 
clear evidence that monocots were both diverse and 
widespread in the middle of the Late Cretaceous 
(Herendeen & Crane, 1995). The dramatic pollen 
character transformations in the Albian-Turonian 
seen in our analysis occurred contemporaneously 
with the rise of angiosperm-dominated ecosystems 
(Lidgard & Crane, 1990; Coiffard et al., 2006). 

The functions of exine architecture in angiosperms 
are considered to be mainly protection (against 
adverse atmospheric conditions such as desiccation 
and UV radiation), storage of physiologically active 
substances, pollen grain clustering, and harmomega- 
thy (Muller, 1979). It 1s likely in the near future that 
we will gain a fuller understanding of Heslop- 
Harrison's (1969, 1972) insights that pollen wall 
patterning is not directly coded by genes but involves 
an interplay between the products of gene expression 
and the influence of physical forces (Blackmore et al., 
2007). Water content status influences not only the 
survival, longevity, and germination ability of 
angiosperm pollen (Franchi et al., 2002; Aylor, 
2003), but also their exine stratification: exine 
reduction or loss in monocots has been associated 
with moist or aquatic environments (Heslop-Harrison, 
1976; Philbrick, 1991; Furness & Rudall, 1999; 
Tanaka et al., 2004), such as those frequented by 
many alismatids (Thanikaimoni, 1978). However, the 
complexity of the exine and the development of its 
stratification under physical stress may limit our 
understanding of the associations between exine 
architecture and habitat moisture. We confirmed that 
three exine stratification traits (exine, tectum, and 
infratectum) are correlated with the evolution of helo— 
hydrophytic habitat (Table 3, Fig. 2). The presence or 
absence of exine is most strongly correlated with 
helo—hydrophytic habitat rather than hydrophytic 
habitat among six exine structures (Table 3, Fig. 2). 
Furthermore, tests of contingent evolution indicated 
that exine absence was not preferentially lost in helo— 
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Appendix 1. 


Taxonomic 


classification 


Stemonaceae 


Velloziaceae 


LILIALES 


Alstroemeriaceae 


Campynemataceae 


Colchicaceae 


Liliaceae 


Melanthiaceae 


Petermanniaceae 
Philesiaceae 


Rhipogonaceae 
Smilacaceae 


ASPARAGALES 
Amaryllidaceae 


Continued. 


Total # of genera 


per family (total # 
of genera sampled) 


16 (2) 


Taxon name, 


with species voucher 


Croomia Torr. 

Stemona Lour. 

S. tuberosa Lour. (Fig. 4C, D) 
Acanthochlamys P. C. Kao 
Talbotia Balf. 


Alstroemeria L. 
Luzuriaga Ruiz & Pav. 
Campynema. Labill. 


Schelhammera R. Br. 
Uvularia L. 


Calochortus Pursh 
Lilium L. 


Tulipa iliensis Regel (Fig. 46— 


I) 
Trillium L. 


T. maculatum Raf. (Fig. 4E, 
F) 

Veratrum L. 

Petermannia F. Muell. 

Philesia Comm. ex Juss. 


Rhipogonum J. R. Forst. & G. 
Forst. 
Smilax L. 


Allium L. 


A. cyaneum Regel (Fig. 4J—L) 


References and voucher information 


References used for palynological data 


Van der Ham, 1991 
Van der Ham, 1991; Furness & Rudall, 2006 


References used for ecological data 


Rogers, 1982; Kubitzki, 1998a 
Rogers, 1982; Kubitzki, 1998a 


China, Yunnan, S. D. Zhang et al. APE-LY-017 (KUN, living collection from KBG) 


Gao, 1987 
Furness € Rudall, 2006 


Sarwar et al., 2010 

Schulze, 1982a; Kosenko, 1994a 

Erdtman, 1952; Dutt, 1970; Goldblatt, 1986; 
Kosenko, 1987 

Kosenko, 1988 

Radulescu, 1973; Buchner & Weber, 2000 
onwards 

Heusser, 1971; Radulescu, 1973; Schulze, 1980b; 
Halbritter & Hesse, 1993; Kosenko, 1999 

Radulescu, 1973; Schulze, 1980b; Halbritter & 
Hesse, 1993; Kosenko, 1999 

China, Xinjiang, D. Y. Tan DY0180 (KUN) 


Takahashi, 1982, 1983; Buchner & Weber, 2000 
onwards 


U.S.A., Oklahoma, R. K. Godfrey 84058 (IBSC) 


Buchner & Weber, 2000 onwards 

Erdtman, 1952; Schulze, 1982a 

Erdtman, 1952; Cranwell, 1953; Heusser, 1971; 
Schulze, 1982a; Cauneau-Pigot, 1988 

Schulze, 1982b; Furness & Rudall, 1999; Chen et 
al., 2006 

Chen et al., 2006; Kodela, 2006 


Nair & Sharma, 1965; Radulescu, 1973; Schulze, 
1980a; Buchner & Weber, 2000 onwards; 
Guler & Pehlivan, 2006 


China, Qinghai, Y. H. Wu LJQ-QLS-2008-0158 (KUN) 


Kubitzki, 1998a 
Kubitzki, 1998a 


Kubitzki, 1998a 
Kubitzki, 1998a 
Kubitzki, 1998a 


Kubitzki, 1998a 
Kubitzki, 1998a 


Kubitzki, 1998a; Dilley, 1999 


Kubitzki, 1998a 


Kubitzki, 1998a; Irwin, 2000 


Mulligan & Munro, 1987; Kubitzki, 1998a 
Kubitzki, 1998a 

Kubitzki, 1998a 

Macmillan, 1972; Kubitzki, 1998a 

Kevan et al., 1991; Kubitzki, 1998a 


Kubitzki, 1998a 
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Appendix 1. Continued. 


Taxonomic 


classification 


Asparagaceae 


Asteliaceae 


Blandfordiaceae 
Boryaceae 


Doryanthaceae 
Hypoxidaceae 
Iridaceae 
Ixioliriaceae 
Lanariaceae 


Orchidaceae 


Tecophilaeaceae 
Xanthorrhoeaceae 


Total # of genera 
per family (total # 
of genera sampled) 


153 (11) 


Taxon name, 


with species voucher 


Agapanthus L'Her. 


Clivia Lindl. 
Asparagus L. 


Arthropodium R. Br. 
Convallaria L. 


Aphyllanthes L. 


Brodiaea Sm. 


Scilla L. 


Anemarrhena Bunge 
Agave L. 
Chlorophytum Ker Gawl. 


Behnia Didr. 
Herreria Ruiz & Pav. 


Astelia Banks & Sol. ex R. Br. 


Blandfordia Sm. 
Alania Endl. 


Borya Labill. 

Doryanthes Corrêa 
Hypoxis L. 

Sisyrinchium L. 
Ixiolirion Fisch. ex Herb. 
Lanaria Aiton 

Neuwiedia Blume 
Cypripedium L. 

Epipactis Zinn 
Tecophilaea Bertero ex Colla 
Xanthorrhoea Sm. 


References and voucher information 


References used for palynological data 


Nair € Sharma, 1965; Radulescu, 1973; Buchner 
& Weber, 2000 onwards 

Afzelius, 1955; Gullvág, 1964 

Nair & Sharma, 1965; Buchner & Weber, 2000 
onwards; Ozler & Pehlivan, 2008 

Radulescu, 1973 

Ma & Hong, 1990 

Radulescu, 1973; Furness, 1985; Lifante et al., 
1990 

Schulze, 1980a 

Radulescu, 1973; Schulze, 1980b; Ghavami et al., 
2009 

Chen & Qiao, 2000 

Alvarez & Kóhler, 1987; Stroo, 2000 

Nair & Sharma, 1965; Buchner & Weber, 2000 
onwards 

Schulze, 1982a; Conran, 1999 

Kosenko, 2001 

Cranwell, 1953; Radulescu, 1973; Kosenko, 2001; 
APSA, 2007 

Schulze, 1982b; Kosenko, 1994b 

Erdtman, 1952; Schulze, 1982b; APSA, 2007 


Erdtman, 1952; APSA, 2007 
Schulze, 1982b; Buchner & Weber, 2000 onwards 
Cranwell, 1953 

Schulze, 1971; Nadot et al., 2006 
Schulze, 1983; Donmez & Isik, 2008 
Erdtman, 1952 

Schill, 1978 

Burns-Balogh & Hesse, 1988 
Ackerman & Williams, 1980 
Erdtman, 1952; Simpson, 1985b 
Chanda & Ghosh, 1976 


References used for ecological data 


Leighton, 1965; Kubitzki, 1998a 


Kubitzki, 1998a; Meerow & Clayton, 2004 
Kubitzki, 1998a 


Kubitzki, 1998a 
Kubitzki, 1998a; Araki et al., 2005 
Kubitzki, 1998 


Han et al., 1991 
Kubitzki, 1998a 


Kubitzki, 1998a 
Gentry, 1982; Kubitzki, 1998a 
Kubitzki, 1998a; Poulsen & Nordal, 2005 


Kubitzki, 1998a 
Kubitzki, 1998a 
Moore, 1966; Kubitzki, 1998a 


Zimmerman & Pyke, 1988; Kubitzki, 1998a 

Kubitzki, 1998a; Western Australian Herbarium, 
1998 onwards 

Keighery, 1984; Kubitzki, 1998a 

Newman, 1929; Kubitzki, 1998a 

Cook, 1990; Kubitzki, 1998a 

Kubitzki, 1998a; Freitas & Sazima, 2003 

Kubitzki, 1998a; Li & Zhou, 2006 

Kubitzki, 1998a 

Benzing & Atwood, 1984; Okada et al., 1996 

Benzing & Atwood, 1984; Davis, 1986 

Brantjes, 1981; Benzing & Atwood, 1984 

Kubitzki, 1998a 

Keighery, 1980; Kubitzki, 1998a 
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Appendix 1. Continued. 


Total # of genera References and voucher information 


Cla 


Taxonomic per family (total + Taxon name, 
classification of genera sampled) with species voucher References used for palynological data References used for ecological data 
Flagellariaceae 1 (1) Flagellaria L. Chanda, 1966; Chanda & Rowley, 1967; Ladd, Linder, 1987; Kubitzki, 1998b 
1977; Linder & Ferguson, 1985 
Joinvilleaceae 1 (1) Joinvillea Gaudich. ex Brongn. Chanda, 1966; Chanda & Rowley, 1967 Kubitzki, 1998b 
& Gris 
Juncaceae 7 (2) Juncus L. Cranwell, 1953; Perveen, 1999: Buchner & Cook, 1990; Kubitzki, 1998b 
Weber, 2000 onwards 
J. sikkimensis Hook. f. (Fig. China, Yunnan, X. X. Zhu APE-LY-037 (KUN) 
6G, H) 
Luzula DC. Cranwell, 1953; Buchner & Weber, 2000 onwards Cook, 1990; Kubitzki, 1998b 
Mayacaceae 1 (1) Mayaca Aubl. Zavada, 1983 Cook, 1990; Kubitzki, 1998b 
Poaceae 668 (2) Anomochloa Brongn. Page, 1978; Michelangeli et al., 2003 Clayton & Renvoize, 1986; Judziewicz & 
Soderstrom, 1989; Givnish et al., 2010 
Oryza, L. Sharma, 1968; Page, 1978; Chaturvedi et al., Clayton € Renvoize, 1986; Givnish et al., 2010 
1998 
Dendrocalamopsis beecheyana China, Yunnan, S. D. Zhang et al. APE-LY-22 (KUN, living collection from KBG) 
(Munro) Keng f. var. 
pubescens (P. F. Li) Keng f. 
(Fig. 6K, L) 
Rapateaceae 16 (1) Rapatea Aubl. Carlquist, 1961; Zavada, 1983; Tiemann, 1985 Kubitzki, 1998b 
Restionaceae 58 (2) Elegia L. Chanda, 1966 Kubitzki, 1998b 
Baloskion Raf. Chanda, 1966; Briggs et al., 1999 Kubitzki, 1998b 
Thurniaceae 2 (2) Prionium E. Mey. Munro & Linder, 1997 Cook, 1990; Kubitzki, 1998b 
Thurnia Hook. f. Zavada, 1983 Cook, 1990; Kubitzki, 1998b 
Typhaceae 2 (2) Sparganium L. Nilsson et al., 1977; Zhang & Chen, 1984 Cook, 1990; Kubitzki, 1998b 
Typha L. Nilsson et al., 1977; Zhang & Chen, 1984; Cook, 1990; Kubitzki, 1998b 
Perveen, 1999 
Xyridaceae 5 (2) Abolboda Bonpl. Carlquist, 1960; Straka & Friedrich, 1984; Cook, 1990; Kubitzki, 1998b; Oriani & Scatena, 
Furness & Rudall, 1999; Campbell, 2012 2011 
Xyris L. Straka & Friedrich, 1984; Rudall & Sajo, 1999; Cook, 1990; Kubitzki, 1998b 
Campbell, 2012 
Outgroups 
AMBORELLALES 
Amborellaceae 1 (1) Amborella. Baill. Sampson, 1993, 2000; Hesse, 2001 Kubitzki et al., 1993; Thien et al., 2003 
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Appendix 2. Comprehensive data matrix of pollen morphological and ecological characters analyzed in this study (see text and 
Table 1 for details). Taxon names are the same as those in Chase et al. (2006) throughout. One hundred twenty monocots and 16 


ç 99 


outgroups are scored, and outgroup taxa are represented in italics. Inapplicable states are coded as “-” and unknown states as “?.” 


Polymorphic states are enclosed in curly brackets. 


Abolboda 00111020---00-0---022 

Acanthochlamys 0100-2010?00211222221 

Acorus calamus 0100-1010100220-10122 

Acorus gramineus 0100-1010100220-10122 

Agapanthus 0100-{1,2}210000210-20?21 

Agave 0100-1211,210(0,1,2100210-20?12,410 

Alania 0100-1110200211120221 

Aletris 0100-11101002{0,1,2}0-{1,2}0121 

Alisma 00211010, 1143 110291110111. 919 

Allium 0100-11,2111010, 1100220-11, 3430120 

Alpinia 09011400===0221 100022 

Alstroemeria 0100-1210:192100220-11, 3,2100 2(0, 122) 
Amborella DIO XOU 11110010, 1H042002 105240021 117014 1, 232 
Anarthria 011110 110110211221211.3]011(1,9) 
Anemarrhena 0100-2110?00210-20220 

Anigozanthos 0201-2212 ,3}00101-0---050 

Anomochloa 011110?10211221000?21 

Aphelia EE 11010, 1Mi0111921031911 

Aphyllanthes 0??11021-13022111??20 

Arisaema 001110{0,1}0---022110-022 

Arthropodium 0100-{1,2}110000210-2??21 

Asarum 0211{0,1}0{1,2¥0,2,3,4}21{0,1}02{1,2}1{0,1,2¥1,2}0{0,1}22 
Ascarina 011001{0,1}10100211020112 

Asparagus 0100-{1,2}110{0,1}002{1,2}0-{1,2,3}012{0,1} 
Asphodelus 0100-1210{0,1}002{1,2}0-{0,1,2}0ł121 

Astelia 0100-{1,2¥0,1}10?002?11?0?2{1,2} 
Austrobaileya 0100-11101002{1,2}0-{2 ,3}0122 

Baloskion 011110?1011122??10112 

Behnia 0100-2210?00210-222211;,2) 

Blandfordia 0100-1110000220-3??5{1,2} 

Borya 0100-{1,2¥0,1}10?00210-2??20 

Brodiaea 0100-2210200210-22221 

Butomus 0100-1110000210-20123 

Cabomba 0100=241,2)1011,2)00221210123 

Calamus 0201-110,1)122000211,2)1(0, 1101, 241,2; 5]0122 
Calectasia 0100-[2,3)1,2)10200210-(1,2]0101 

Calochortus 01009211,211,210000210-—2292[05 1,2] 
Calycanthus 0200-10. 110092101, 0, 2) 100211, 230-11, 935 510122 
Campynema 0100-2110?00210-{1,2}0??1 

Canna 00111020- --00-0---04{4,5}2 

Carex EE e 2 624 0652119020 ee EH 
Carludovica 01112{0,1}111?1?220-10122 

Cartonema 0110-1210100210-6??21 

Ceratophyllum 00111010- --00-0-- -033 

Chloranthus 0200-010,1)42110,11021(0, 11120122 
Chlorophytum 0100-[1,21110100220-1??[0,2]0 

Chorigyne 0100-??107?00220-1?122 

Clivia 0100-2210000210-20?{2,5}0 

Convallaria 0100-{1,2}110000220-10122 

Costus 0?211{0,1}3{1,4}30{1,3}01-0-0-0{2,5}2 
Croomia 0100-1110100210-20121 

Cyclanthus OTITIO. D ITO, 11 1001022010122 

Cymodocea 00-2--40---00-0-- -033 

Cypripedium 0100-1110200220-{0,1}0021 

Dasypogon 0100-1110000220-101{2,5}1 


Dioscorea 011,2)01-11,2,3H0,1H1,2,3)21002(1,2)0-(1,2,3,4J)012(1,2) 
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Doryanthes 
Drimys 
Ecdeiocolea 
Elegia 
Epipactis 
Eriocaulon 
Euterpe 
Flagellaria 
Freycinetia 
Gymnostachys 
Hanguana 
Heliconia 
Hemerocallis 
Herreria 
Hypoxis 
Illicium 
Ixiolirion 
Japonolirion 
Joinvillea 
Juncus 
Kingia 
Lactoris 
Lanaria 
Lilium 
Liriodendron 
Luzula 
Luzuriaga 
Magnolia 
Mapania 
Maranta 
Mayaca 
Murdannia 
Musa 
Narthecium 
Neuwiedia 
Nymphaea 


Nypa 
Orchidantha 
Orontium 
Oryza 
Pandanus 
Petermannia 
Petrosavia 
Philesia 
Philydrum 
Pleea 
Pontederia 
Potamogeton 
Prionium 
Puya 
Rapatea 
Rhipogonum 
Sagittaria 
Saururus 
Schelhammera 
Schisandra 
Scilla 
Sisyrinchium 
Smilax 


Sparganium 


Continued. 
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0100-{1, 2}110200210-2727{2,5}1 
111100{0,1}10111210-2012{1, 2} 
011170110211221210111 
011110110?11220-1?112 
11111011010,1)10210-2?122 
05211010,114/1.2,91213022110,110,1,3]1012]1,2] 
0100-2{1,2}101002{1,2}0-{1, 2}0?2{1, 2} 

Old 11 CIAO LI 25 11:10 11015927 
0111{0,1}1{0,1}1001122{0,1}0{0,1}01{4,6}2 
0100-1110100220-{1,5}0122 
001110{0,1}0---00-0---022 
0111{0,1}020---00-0---0{4,5}2 

0100-{2, 3}{2,3}10000210-20121 
0100-?110?00210-207??2 

0100-{1, 2}110100210-20?2{1, 2} 
0211{0,1}01321{0,4}0210-20121 
0100-211,2)10200210-20?20 
0100-1010?00200--01?2 
Og111011011122]7011011 

TIi11010. 1140110223070 ,1]921]] 9101.2] 
0100-1{1,2}10000220-{0,1}01{2,5}1 
11{0,1}{0,1}{0, 1}0{0, 1}101{0,1}0210-{1,3}{0,1}111 
0100-1{0, 1}10000220-10121 
O1{0,1}0-{1,2}{2,3}{1,2,3}01{0,1}0210-20121 
0100-{1,2}{1,2}101002{1,2}0-{1,2,3}012{1,2} 
LTC LOTTO E 2,27 
0100-1{0, 1}10100210-27702 
0100-1246, 2101002210 1300 D 1, 2H 241, 2} 
111110101) 101 T1211; 23404 10 19097 SE 
00111020-- -00-0- - -022 
0100-1?10100210,1]0-20122 
0100-211010022110,2)10121 
001110(2,3]0- - -00-0-- -014,5)2 

0100-1{0, 1}102002{1,2}0-{1,2}0122 
0100-1010200210-20122 

OFT 210 209, 1340, IO, 1H 1270; 1140 2022110; 2110, 2110; 1 5 122 
021110112020211,2]11111,2]0122 
01012030-- -00-0- - -022 
0100-{1,2}210100220-{1,5}0122 
011110110211211{0,1}{0,6}0?711 
01110:052)01052)1200110& 1171621414 DO MO 22040714 2,507210; 29022,4212 
0100-2110?0022115??22 
0100-21[0,1)10100210-11,210?(0,2M1,2] 
00110020- - -021110??511,2] 
1100-1110000210-20122 
0201-1[0,1)2??00210-20121 
0100-21{2,3}{0,2}000220-{0,3,5}0123 
{0,1}011{0,1}0{0, 1}0---0211120113 
11111001011022100??12 
0100-1210000210-201(4,5)1 

022 1101121202101,210-11,2]912(1.,2) 
0100-(1,2)110?0021(0,1)120122 

002110{0, 1}343110221{0,1}101{1,2}3 
0100-2010100221010122 
0100-11,2)210?002????0?22 
011100{0,1}{3,4}{0, 2}1{0,4}0210-20121 
0100-{1,2}1100002{1,2}0-{1,2}0121 
0100-2{1,2}10000210-20121 

001110{0, 1}0---02{0, 2}1131022 
011110110110210-20?12 
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Appendix 2. Continued. 


Sphaeradenia 0110,140,1H0, 112110010, 110211, 2}0-{1, 2}012{1, 2} 
Stemona 0100-1{0, 13101002{0,1,2}0-{0,1,2,3,5}{0,1,2}12{0,1} 
Strelitzia 00111030---00-0---052 

Tacea 0100-{1, 2}1101002{1,2}0-{2,3}01{0, 2}2 

Talbotia 0100-1010200220-301?2 

Tecophilaea 0100-1110200220-301?1 

Thurnia 111110010110211160?12 

Tillandsia O100-{1, 2}{1,2}1{0,1}{1,2}002{1, 2}0-{1,2}01{2,5}{0,1} 
Tofieldia 0201-310; 11221002110, 1170, 21112101211, 21 

Trichopus 0100-10101002211???02 

Triglochin 00111010, 110- - -0210-20112 

Trillium 00111010, 110- - -0(0,2H0, 110-3010, 1}2{1, 2} 

Typha 10131 131010; 1310110210=20212 

Uvularia 0100-(1,2+1(1,2)0200211(0,1)10?2(1,2) 

Vallisneria 00111030---00-0---033 

Veratrum 0100-11,2)110200210-2??22 

Xanthorrhoea 0100-11,2)110?00210-201(2,5)1 

Xeronema 0100-1110000210-2?1(2,5y0 

Xyris 0100-11{1,2}0{1,2}00210-2772{1,2} 


Zostera 00-2--40---00-0---033 
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Democratic data matrix of pollen morpholog- 


ical and ecological characters analyzed in this study (see text 


and Table 1 for details). Taxon names are the same as those in 
Chase et al. (2006) throughout. One hundred twenty monocots 
and 16 outgroups are scored, and outgroup taxa are represented 


in italics. Inapplicable states are coded as 


states as “?.” 


Abolboda 


Acanthochlamys 
Acorus_calamus 
Acorus_gramineus 


Agapanthus 
Agave 

Alania 
Aletris 
Alisma 
Allium 
Alpinia 
Alstroemeria 
Amborella 
Anarthria 
Anemarrhena 
Anigozanthos 
Anomochloa 
Aphelia 
Aphyllanthes 
Arisaema 
Arthropodium 
Asarum 
Ascarina 
Asparagus 
Asphodelus 
Astelia 
Austrobaileya 
Baloskion 
Behnia 
Blandfordia 
Borya 
Brodiaea 
Butomus 
Cabomba 
Calamus 
Calectasia 
Calochortus 
Calycanthus 
Campynema 
Canna 

Carex 
Carludovica 
Cartonema 
Ceratophyllum 
Chloranthus 
Chlorophytum 
Chorigyne 
Clivia 
Convallaria 
Costus 
Croomia 
Cyclanthus 
Cymodocea 
Cypripedium 


ee 99 


and unknown 


0O111020>=00=0===0 
O100=2010 200211227? 
O100=1010100220-101 
0100-1010100220-101 
0100—2210000210-207 
0100-1210100210-20? 
Q100-1110?00211120? 
0100-1110100210-201 
0021100431102211101 
0100-1110100220-401 
001110.20>==022 11000 
0100-1210100220-30? 
0111000102102111701 
0111101102112212301 
0100-2110?00210-20? 
0200-22200101-0---0 
011110?10211221000? 
OIl1l10010111221031? 
O 22 11021=13022141122 
00111000---022110-0 
0100-2110000210-2?? 
0211102321002112201 
0110010101002110201 
0100-2110100220-101 
0100-1210100220-101 
0100-1010?002?11?0? 
0100-1110100210-201 
011110?1011122??101 
0100-?210?00210-2?? 
0100-1110000220-3?? 
0100-2110?00210-2?? 
0100-2210?00210-2?? 
0100-1110000210-201 
0100-22102002212101 
0201-1122000210-201 
0100-2210200210-201 
0100-2110000210-2?? 
0200-0122100210-201 
0100-2110?00210-20? 
00111020---00-0---0 
1121101401102210121 
011120111?1?220-101 
0110-1210100210-6?? 
OOTT11O10==>=00=0===0 
0200-00421002101201 
0100-2110100220-1?? 
0100-27210700220-1?1 
0100-2210000210-20? 
0100-2110000220-101 
0?21103430101-0-0-0 
0100-1110100210-201 
011101110010220-101 
00-2--40---00-0---0 
0100-1110200220-000 
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Appendix 3. 


Dasypogon 
Dioscorea 
Doryanthes 
Drimys 
Ecdeiocolea 
Elegia 
Epipactis 
Eriocaulon 
Euterpe 
Flagellaria 
Freycinetia 
Gymnostachys 
Hanguana 
Heliconia 
Hemerocallis 
Herreria 
Hypoxis 
Illicium 
Ixiolirion 
Japonolirion 
Joinvillea 
Juncus 
Kingia 
Lactoris 
Lanaria 
Lilium 
Liriodendron 
Luzula 
Luzuriaga 
Magnolia 
Mapania 
Maranta 
Mayaca 
Murdannia 
Musa 
Narthecium 
Neuwiedia 
Nymphaea 
Nypa 
Orchidantha 
Orontium 
Oryza 
Pandanus 
Petermannia 
Petrosavia 
Philesia 
Philydrum 
Pleea 
Pontederia 
Potamogeton 
Prionium 
Puya 
Rapatea 
Rhipogonum 
Sagittaria 
Saururus 
Schelhammera 
Schisandra 
Scilla 


Sisyrinchium 


Continued. 


0100-1110000220-101 
0200-1122100220-101 
0100-2110200210-2?? 
111100110111210-201 
0111101102112212101 
011110110?11220-1?1 
111110110010210-2?1 
0??11011?1302211101 
0100-2110100220-10? 
0111101101112212301 
011101010011220-101 
0100-1110100220-101 
00111000- --00-0---0 
0O11110205-5-00-0-&-0 
0100-3310000210-201 
0100-?110?00210-20? 
0100-2110100210-20? 
021110132100210-201 
0100-2110200210-20? 
0100-1010?00200--01 
0111101101112210110 
IIITIOQIOITO2210477 
0100-1210000220-101 
111110110100210-311 
0100-1110000220-101 
0100-2210100210-201 
0100-2210100220-101 
11111011011022101?? 
0100-1110100210-2?? 
0100-2210100220-101 
1111100101112210101 
001110?0---00-0---0 
0100-1?101002-0-201 
0100-21101002212101 
00111020- --00-0-- -0 
0100-1010200210-201 
0100-1010200210-201 
0211101120202212011 
0211101120202211101 
01012030- --00-0---0 
0100-2210100220-101 
011110110211211060? 
0111200100101-0-- -0 
0100-2110?0022115?? 
0100-2010100210-20? 
00110020---0211102? 
1100-1110000210-201 
0200-102??00210-201 
0100-2122000220-301 
00111000---02111201 
11111001011022100?? 
0100-1210000210-201 
021110112120210-201 
0100-1110?00210-201 
0021100431102211101 
0100-20101002210101 
OLOO=2210 72002777207 
011100142140210-201 
0100-2110000210-201 
0100-2110000210-201 
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Smilax 
Sparganium 
Sphaeradenia 
Stemona 
Strelitzia 
Tacca 
Talbotia 
Tecophilaea 
Thurnia 
Tillandsia 
Tofieldia 
Trichopus 
Triglochin 
Trillium 
Typha 
Uvularia 
Vallisneria 
Veratrum 
Xanthorrhoea 
Xeronema 
Xyris 
Zostera 


00111000---0200--10 


011110110110210- 20 2 
0101-21 100004220-1091 
0100-1110100200--11 


OOIITOUSO0---00-0-- =O 


0100-2110100220-301 
0100-1010200220-301 
0100-1110200220-301 
111110010110211160? 
0100-1110100210-201 
0200-1022100210-201 
0100-10101002211??? 


00111000---0210-201 
001110105 ==00=-0===0 


111110110110210-20? 
0100-1110200211110? 


00111030---00-0---0 


0100-2110200210-2?? 
0100-1110?00210-201 
0100-1110000210-2?1 
0100-1110200210-2?? 


00-2--40---00-0---0 


Annals of the 


Missouri Botanical Garden 


Appendix 4. Data matrix of pollen morphological and 


ecological characters used in analyses of correlated evolution 


(see text and Table 1 for details). Taxon names used are the 
same as those in Chase et al. (2006) throughout. One hundred 
twenty monocots and 16 outgroups are scored, and outgroup 


taxa are represented in italics. Unknown states are coded as 


ce) KM 


Abolboda 
Acanthochlamys 
Acorus_calamus 
Acorus_gramineus 
Agapanthus 
Agave 

Alania 

Aletris 
Alisma 
Allium 
Alpinia 
Alstroemeria 
Amborella 
Anarthria 
Anemarrhena 
Anigozanthos 
Anomochloa 
Aphelia 
Aphyllanthes 
Arisaema 
Arthropodium 
Asarum 
Ascarina 
Asparagus 
Asphodelus 
Astelia 
Austrobaileya 
Baloskion 
Behnia 
Blandfordia 
Borya 
Brodiaea 
Butomus 
Cabomba 
Calamus 
Calectasia 
Calochortus 
Calycanthus 
Campynema 
Canna 

Carex 
Carludovica 
Cartonema 
Ceratophyllum 
Chloranthus 
Chlorophytum 
Chorigyne 
Clivia 
Convallaria 
Costus 
Croomia 
Cyclanthus 
Cymodocea 


00000010 
11011000 
11111011 
11111011 
11011000 
11011000 
11011000 
11011000 
11111011 
11111000 
11111000 
11111000 
11111000 
11111000 
11011000 
10000000 
11111000 
11110100 
11111000 
11111000 
11011000 
11011000 
11011000 
11011000 
11011000 
11?11000 
11011000 
11111010 
11011000 
11111000 
11011000 
11011000 
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